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ABSTRACT 


We  study  the  initial  value  problem  for  the  nonlinear  Volterra  integro- 
differential  equation 


(+) 


-py 


' (t ) = J a(t  - s)F(y(t),y(s))ds  (t  > 0) 


y(t)  = g(t) 


(-00  < t < 0)  , 


where  p > 0 is  a small  parameter,  a is  a given  real  kernel,  and  F,  g 
are  given  real  functions;  (+)  models  the  elongation  ratio  of  a homogeneous 
filament  of  a certain  polyethylene  which  is  stretched  on  the  time  interval 
(--w,0]  , then  released  and  allowed  to  undergo  elastic  recovery  for  t > 0. 

Under  assumptions  which  include  physically  interesting  cases  of  the  given 
functions  a,  F,  g,  we  discuss  qualitative  properties  of  the  solution  of  (+) 
and  of  the  corresponding  reduced  problem  when  p = 0,  and  the  relation 
between  them  as  p 0 , both  for  t near  zero  (where  a boundary  layer  occurs) 
and  for  large  t.  In  particular,  we  show  that  in  general  the  filament  does  not 
recover  its  original  length,  and  that  the  Newtonian  term  -py'  in  (+)  has 
little  effect  on  the  ultimate  recovery  but  significant  effect  during  the  early 
part  of  the  recovery. 
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1.  Introduction.  We  study  the  nonlinear  Volterra  integrodifferential  equation 

t 

(i.  1)  -py'(t)  = / a(t  - s)F(y(t),y(s))ds  (t  > 0;  ' = d/dt) 

_ 00 


subject  to  the  initial  condition 

(1.2)  y (t ) = g(t)  (-°°  < t < 0 > . 

This  initial  value  problem  arises  as  a mathematical  model  for  a process 
in  polymer  rheology  which  is  described  in  Appendix  A.  In  the  specific 
problem  discussed  there  p is  a positive  parameter  related  to  part  of  the 
viscosity,  and  the  given  real  functions  a,  F,  g take  the  forms 

m 


(1.3) 


a(t)  = Tj  a exp(-t/ t ) 
k = l k 


where  a^  and  are  positive  constants, 


(1.4) 

and 


3 2 

F(y,z)  = y /z  - z , 


1 if  -°°  < t < - tQ 


(t0  > 0) 


(1.5) 


g(t)  = 


K(tn  ) 

e if  -tQ  < t < 0 , 
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where  k is  a positive  constant.  The  unknown  function  y measures 
the  ratio  of  the  extended  length  to  the  original  length  of  a homogeneous 
filament  which  is  stretched  in  accordance  with  (1.2)  on  the  time  interval 
t-30,!)]  and  then  released.  The  equation  (1.  1)  then  describes  the  process 
of  elastic  recovery.  Also  of  interest,  both  mathematically  and  physically, 

I 

is  the  reduced  equation 

t 

(1.6)  0 “ / a<t  ' s)F(y(t),y(s))ds  (t  > 0)  , 

-00 

with  y(t)  = g(t)  on  (-°°,0),  and,  in  particular,  the  relation  between 
the  solutions  of  ( 1 . 1 ) , (1.2)  and  of  ( 1 . 6 ) as  p — 0+  ; for  small  p > 0, 

(1.1),  (1.2)  may  be  regarded  as  a singular  perturbation  of  (1.6). 

The  purpose  of  this  paper  is  to  discuss  the  qualitative  behaviour  of 
solutions  of  ( 1 . 1),  (1 . 2)  on  the  one  hand  and  of  (1.6)  on  the  other,  and 
the  relation  between  them  as  p — 0+ . Our  analysis  is  not  confined  to 
the  specific  forms  of  a,  F,  g in  (1.3)  -(1.5),  but  rather  we  abstract 
the  essential  properties  of  these  functions.  (See  also  (A28)  - (A33)). 

Two  results  of  particular  interest  are:  (i)  in  general,  the  filament  does 
not  return  to  its  original  length,  as  confirmed  by  experiments  (see  Appendix  A 
and  Theorems  3 and  4);  (ii)  similarly  to  behaviour  in  singular  perturbation 
problems  for  ordinary  differential  equations,  the  solution  of  (1.1),  (1.2)  for 
small  p > 0 decreases  rapidly  as  p decreases  to  the  solution  of  (1.6) 
near  t = 0,  becoming  close  to  it  in  a "boundary  layer"  time  interval  of 
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order  j±  |log  (jl  I , and  thereafter  remains  close  (Corollary  2.  1 and 
Theorems  7 and  8).  The  introduction  of  the  Newtonian  term  -py'  in 
(1.1)  has  little  effect  on  the  ultimate  behaviour  of  the  solution. 

From  the  mathematical  point  of  view  the  theory  of  integrodifferential 
equations  with  decreasing  convex  kernels  and  with  nonlinearities 
consisting  of  functions  of  one  variable  is  well  known,  see  e.g.  [7]  , [18] 
where  references  to  other  literature  are  given.  A part  of  the  novelty  of 
the  present  analysis  is  that  F in  ( 1 . 1 ) and  (1.6)  is  a function  of  two 
variables,  y(t)  and  y(s). 

We  acknowledge  with  pleasure  the  various  helpful  discussions  with 
colleagues  during  the  preparation  of  this  paper,  in  particular  with 
M.  G.  Crandall,  R.  W.  Dickey,  F.  Howes,  J.  J.  Levin,  S.  O.  Londen, 

S.  V.  Parter,  A.  Pazy,  D.  F.  Shea,  O.  J.  Staffans,  L.  Tartar,  and 
W.  Wasow. 
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2.  Statement  of  Results.  Let  IR  denote  the  real  numbers,  IR  the 
positive  real  numbers,  and  C the  set  of  k times  continuously 
differentiable  functions. 

We  make  the  following  assumptions  on  the  functions  a,  F,  g 


throughout: 


a*  C [0,'»);a(t)>0,  a 1 (t>  < 0 (0  < t < °o)  : 
a«  L1  (0, »);  log  aft)  is  convex  (0<t<°°),  i.e. 
a'(t)/a(t)  is  nondecreasing: 


F : IR  x IR  — IR;  F(x,x)  - 0 for  every  x > 0 ; 

F e cNlR  X IR  ) and  F^(y,z)  > 0,  F^(y,z)  < 0 

(y,z  t IR+),  where  the  subscripts  denote  partial  differentiation; 


j g : - IR  ; g(-°°)  = 1,  g(0)  > 1 ; 

(H  ) < 

9 1 g t C(-°°,0]  and  g is  nondecreasing  . 

It  is  readily  verified  that  the  specific  functions  defined  in  (1.3)  - (1.5) 
satisfy  (H  ),  (H  ) and  (H  ) respectively.  While  for  the  majority  of 

3 I-  CJ 

results  these  are  the  only  conditions  required,  some  additional  assumptions 
are  needed  in  Theorems  3 and  4,  Corollary  5.  1,  and  Theorems  6 and  8. 

The  first  result  concerns  the  global  existence  and  uniqueness  of 
the  solution  of  (1.  1),  (1.2)  for  a fixed  p > 0 and  gives  some  useful 
properties  of  the  solution. 

Theorem  1.  Let  (H  ).  (H  ).  (H  ) be  satisfied.  Then  for  each  p > 0,  the 
initial  value  problem  (1.  1).  (1.2)  has  a unique  solution  4>(t,|j.)  on  [0,°°) 
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(2.  1) 


(2.  2) 


fylng  the  following  properties : 

4 ' (t , pt)  < 0 and  1 < <t  (t,  p)  < g(0)  (0<t<oo); 

if  gj,g2  satisfy  (H^ ) and  if  g j (t ) > g2(t)  (-<»  < t < 0), 

\ then  the  corresponding  solutions  4 j (t , pu) , 4 2(t , pt)  of  (1.1), 
(1.2)  satisfy  <fc  j (t,  p)  > 4 ^(t , p)  (0  < t < °o  ) ; 


(2.  3) 


{ 


if  (Xj  > 
satisfy 


P2>  then  the  corresponding  solutions  of  ( 1 . 1 ) . 
4>(t,p  ) > 4>(t, p ) (0  < t < 00  ) . 


(1.2) 


An  immediate  consequence  of  (2.1 ),  (Z.  3)  is 

Corollary  1.1.  Let  (H  ).  (H„),  (H  ) be  satisfied.  Then  for  each  fixed 
a f g 

p > 0 one  has 

(2.4)  a(p)  = lim  4 ( t , p ) exists  and  a(p)  > 1 ; 

t—  °° 

moreover.  If  p^  > p^,  then  a(p^)  > «(p  ) > 1.  Theorem  1 is  proved  in 
Section  3. 

Remark  1.1.  In  the  special  case  a (t)  s Ae  ^ A > 0,  X.  > 0 which 

satisfies  (H  ) (with  a'^  = -\),  conclusions  (2.1)  and  (2.4)  can  be 
a a(t) 

strengthened  respectively  to  the  following: 

(2.1')  4' (t, p)  < 0 and  1 < yQ  < 4(t,  p)  < g(0)  (0  < t < «> ) , 

(2.4')  a-(p)  = lim  4(t , p)  exists  and  q'(pl)  > y > 1 , 

t -►  00 

where  (see  the  proof  of  (2.  5)  in  Theorem  2)  yQ  is  uniquely  defined  by 
the  equation 
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J eXsF(yQ,g(s))ds  = 0 . 

-00 

The  proof  of  Remark  1.1  is  carried  out  in  the  course  of  proving  (2.1)  (Case  (ii)). 
Remark  1.2.  If  one  is  interested  only  in  global  existence,  rather  than 
further  properties  of  solutions  of  (1.1),  (1.  2),  then  the  following  existence 
result  can  readily  be  established. 

Proposition  1.  Let  a c Ll(0, »),  a(t)  > 0 (0  < t < «> );  let  g satisfy  (H  );  jet 
F t C(IR+ x F+),  F(x,x)  = 0 for  every  xe  E+  and  F(y,z)<0  for 

t 

y < z and  F(y , z)  > 0 for  y > z (0  < y , z < ) . Then  for  each  fixed  p > 0 

the  initial  value  problem  (1.1),  (1.2)  has  at  least  one  solution  on  0 < t < °o  . 
The  proof  of  this  result  follows  from  the  following  observations:  (i)  under 
the  present  hypotheses  (1. 1),  (1.  2)  has  a local  solution  on  0 < t < t^ 

(see  e.g.  [17,  Lemma  1.1]  , also  [16]  );  (ii)  if  y(t)  is  any  solution  of  (1.1), 

(1.  2)  on  0 < t < oo  , then  a slight  modification  of  the  early  part  of  the 
proof  of  Theorem  1 shows  that  1 < y(t)  < g(0)  (0  < t < « );  (iii)  in  view 

of  (i),  (ii),  every  such  local  solution  can  be  extended  (but  not  necessarily 
uniquely)  to  the  interval  [0,°°)  (see  e.g.  [17,  Lemma  1.  2 ],  also  [ 16] ) . 

Concerning  the  reduced  problem  (1.6)  we  prove  the  following  global 
results . 

Theorem  2.  Let  (H  ).  (HJ,  (H  ) be  satisfied.  Then  (1.6)  has  a unique 
a r g 

(continuous)  solution  on  [ 0 , °°  ) satisfying  the  following  properties : 

(2.5)  if  yQ  - 4>0<°)>  foen  1 < yQ  < 9(°)  : 


if  aft)  ^ Ae  , A > 0,  X.  > 0,  then 

(2.6)  $0  c C*[  0 , °°  ),  c^(t)  < 0 and  1 < q.  (t)  < yQ  (0  < t < oo 

if  aft)  £ Ae  A > 0,  X > 0,  then  4>Q(t)  = yQ  (0  < t < °°  ); 


(2.7) 


if  gj,  g 2 satisfy  (Hg)  and  if  g^t)  > g2(t)  ( -oo  < t < 0) , 
then  the  corresponding  solutions  4>q^,  4>q2*  of  (1.6) 
satisfy  4^  *(t)  > 4>Q2'(t)  (0  < t < °°  ) ; 


r 

ii  <t>(t , |jl ) is  the  solution  of  (1.1).  (1.2)  for  a fixed  p > 0, 

(2.8)  < and  if  <J>Q(t)  is  the  solution  of  (1.6).  then  $ (t)  < 4>(t,p) 

^ (0  < t < 00  ) . 

As  a consequence  of  (2.4),  (2.6),  (2.8)  we  have  the  first  statement  in 
Corollary  2.1.  Let  (H^).  (H^),  (H  ) be  satisfied.  Then 

(2.9)  a = lim  cj>  (t)  exists  and  1 < a < a(p)  ; 

t—  oo 

if  a (t)  = Ae  *’t,  A > 0,  X.  > 0,  then,  in  fact,  = y > 1.  Moreover. 

(2.10)  lim  ar(p)  = or  . 

Theorem  2, and  (2.10)  are  proved  in  Section  4.  In  Theorem  8 below  we 
establish  a more  precise  result  than  (2.10)  under  some  additional  assumptions. 

The  next  task  is  to  establish  the  physically  important  fact  that  the 
limiting  value  o-(p)  of  the  solution  4>(t,  (jl)  of  (1.1),  (1.  2)  as  t - oo 
satisfies  ar(p)  > 1 (p  > 0),  rather  than  the  weak  form  c*(p)  >1  in  (2.  4). 
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By  properties  (2.8),  (2.9)  it  suffices  to  prove  > 1 (aQ  - lim  <j,  (t)). 

t * to 

Theorem  3.  Let  (H  ) (H  ),  (H  ) be  satisfied.  If.  in  addition. 

a ’ p > g — » > 

00 

(2. 11)  J ta(t)dt  < 00  , 

0 


then 

(2.12)  tto>1' 

We  remark  that  in  the  special  case  a(t)  = Ae  Kt,  A > 0,  X.  > 0,  there  is 
nothing  to  prove  since  ct>0 (t ) sy  >1,  Theorem  3 is  proved  in  Section  5. 

Theorem  3 is  best  possible  in  the  following  sense. 

Theorem  4.  Let  (H^).  (H  ) be  satisfied  and  let 

(2.13)  F(y , z)  = y - z . 

If 


00 

(2.14)  J sa(s)ds  = 00  , 

0 

then 

(2.  1 5)  lim  <t  (t)  = 1 , 

t-  00 

when  d>Q  is  the  solution  of  (1.6)  with  F satisfying  (2.  13). 

Remark  4.1.  F(y,z)  = F^fy  - z),  where  FQ  > 0 is  a constant,  is  the 
simplest  form  of  F being  linear  consistent  with  assumptions  (H  ). 

Remark  4.2.  A similar  proof  applied  to  (1 . 1),  ( 1 . 2)  shows  that  if  the 
hypotheses  of  Theorem  4 hold,  then  lim  4>(t , |x)  = 1,  where  <f(t,p)  is 

t—  oo 

the  solution  of  ( 1 . 1) , (1.2)  with  F satisfying  (2.  1 3) . Theorem  4 is  proved 
in  Section  6 . 
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We  now  turn  to  the  study  of  asymptotic  behaviour  as  t -*•  of 
solutions  of  ( 1 . 1 ) , ( 1 . 2)  and  of  the  reduced  problem  (1.6).  The  first 
result  is  a useful  estimate  for  <t’(t,p),  where  <f  is  the  solution  of 
(1  . 1),  (1.2). 

Theorem  5.  Let  (H  ).  (H„),  (H  ) be  satisfied.  Then  there  exists  a 
a F ’ g 

constant  K(  > 0 (independent  of  p)  and  constants  pQ  > 0 and 
K K(pQ)  > 0 such  that  the  solution  <t(t,p)  of  (1.1).  (1.2)  satisfies 
the  estimate 

~ 00 

K ^ 

(2.  16)  0 < - cF ' (t , jjl)  < — exp(-K.t/p)  + K f a(s)ds  (0<t<£»;0<p<  p ) . 

~ p 1 J ~ 0 

As  a consequence  of  Theorem  5 and  the  logarithmic  convexity  of  a,  one 
obtains 

Corollary  5.1.  If,  in  addition.  (2.11)  is  satisfied,  then  there  exist 
constants  pQ  > 0 and  K = K(p  ) > 0 such  that 

00 

(2.17)  0 < <t(t,p)  - a(p)  < K / (s  - t)a(s)ds  (0  < t < »;  0 < p < p ) , 

t U 

where  ar(p)  = lim  <t(t,p)  > 1 (by.  Theorems  1 and  3).  Moreover, 
t—  oo 

00 

(2.18)  0 < <f(t,p)  - «t»Q(t)  < ^(p)  - aQ  + K f (s  - t )a(s)ds  (0  < t < oo;  0 < p < pQ)  , 

where  is  the  solution  of  (1.6)  having  aQ  = lim  4.  (t)  > 1 (by 

t-*  00 

Theorems  2 and  3). 

Theorem  5 and  Corollary  5.1  are  proved  in  Section  7. 
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Remark  s.l.  An  estimate  similar  to  (2.16)  holds  for  the  solution  A of 

0 

k' 

the  reduced  problem  (1.6),  (but,  of  course,  without  the  term  - exp(-K,t/p)). 

M-  1 

This  can  be  proved  as  a special  case  of  the  proof  of  Theorem  5 by  obtaining 
an  estimate  of  the  form  (7.6)  in  Section  7,  where  f is  now  independent 
of  pi,  and  using  it  and  an  estimate  of  the  form  (7.4)  (also  now  independent 
of  p)  in  (4.2)  of  Section  4.  Consequently,  (2.17)  also  holds  with  <j>(t,p) 
replaced  by  dQ(t)  and  a(p)  by  a and  with  0 < t < eo  . 

In  connection  with  the  last  statement  in  Remark  5.1  it  is  of  interest 
to  note  that  this  estimate  for  cQ  can  be  proved  independently  of  Theorem  5, 
and  we  state  this  fact  as  a separate  result  in  Theorem  6.  However,  it 
should  also  be  noted  that  the  estimate  (2.16)  (for  $' ) described  in 
Remark  5.1  cannot  be  obtained  from  Theorem  6. 

Theorem  6.  LetlH^),  (Hp),  (H  ) be  satisfied,  and  let  (2.11)  hold.  Then 
there  exists  a constant  K > 0 such  that 

00 

(2.19)  0 < dQ(t)  - <*0  < K J (s  - t)a(s)ds  (0  < t < <»)  . 

Theorem  6 is  proved  in  Section  8. 

Remark  6.1.  One  can  also  prove  the  estimate  (2.17)  for  the  solution  d(t,p) 
of  (1.1),  (1.2)  in  the  manner  of  Theorem  6 , without  using  (2.16). 

Our  next  task  is  to  establish  the  existence  of  a boundary  layer  in  a 
neighbourhood  of  t = 0 as  p -»  0 . For  this  purpose  we  consider  the 
following  approximation  of  the  problem  (1.1),  (1.  2)  for  small  t > 0; 

0 

(2.20)  -pv'(t)  = / a(-s)F(v(t) , g ( s ) )ds  (t  > 0;  v(0)  = g(0))  . 

- oo 


L 
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It  will  be  observed  that  (2.  20)  is  not  a Volterra  equation,  but  acts  rather 
like  an  ordinary  differential  equation.  Performing  the  stretching  transformation 
(2.  21)  t = jit 

and  setting  w(t)  v(t)  transforms  (2.20)  to 
, 0 

(2.22)  - — i a(-s)F(w(T),g(s))ds  (t  > 0;  w(0)  = g(0))  . 

dT  .00 

Theorem  7.  Let  (H  ),  (H_),  (H  ) be  satisfied.  Then  the  initial  value 
problem  (2.  22)  has  a unique  solution  w = £(t)  existing  on  0 < t < oo 
and  satisfying  the  following  properties: 

(2.  23)  lim  |(t)  - yQ  = ^(O);  0 < 4 (r)  - yQ  < (g(0)  - yQ)e  Kt  (0  < t < « ) , 

T - » 

where  6^  is.  the  solution  of  (1.6)  (see  Theorem  2)  and  K is  some  positive 
constant. 

Moreover,  if  <f>(t,  p)  is  the  unique  solution  of  (1. 1),  (1.2)  ira 
Theorem  1 and  if  |(t/p)  is  the  unique  solution  of  (2.  20)  for  p>0,  then 
for  anv  t > 0 there  exists  a constant  K > 0 (independent  of  p) 

i-  o 

such  that 

(2.24)  U(t,p)  - |(t/p)|  < Kt+0(p)  (0<  t<  tQ:  p-  0f)  . 

The  estimate  (2.  24)  establishes  the  existence  of  a boundary  layer  in  a 
positive  neighbourhood  of  t = 0.  Theorem  7 is  proved  in  Section  9. 

In  Corollary  2.1  we  showed  that  a(p)  — aQ  as  p -►  0 , so  that 
the  solutions  $(t,p)  of  (1. 1) , (1.  2)  and  <»0(t>  of  (1.  6)  do  not  differ  by 
much  for  small  p > 0 and  for  large  t.  Our  final  result  makes  this 
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precise,  under  the  additional  assumptions  that  ta(t)  e L (0,ao  ) and 
F < C2I1R'  x ]r‘  ). 

Theorem  8.  Let  (H  ).  (H„),  (H  ) be  satisfied.  In  addition,  assume  that 

a } ’ g 

2 4 + 

Ft  C (1R  x 1R  ) and  that  (2.11)  holds.  Then  there  exist  constants  K > 0, 

> 0 and  a function  -y  t C^f  0 , °°  ),  y positive . bounded  and  non- 
decreasing.  such  that 

(2.  25)  cQ(t)  < c( t,p)  < bQ(t)  + (9(0)  - <t>0(0))exp(-Kt/p)  + y(t)p  |log  p | 

(0  < t < oo  ; 0 < |jl  < p ) . 

In  particular,  as  an  immediate  consequence  of  (2.  25) , there  exists  a 
consta nt  K > 0 such  that 

(2.  26)  0 < c(t,p)  - <t Q (t)  = 0(p  |log  pl),(p  - 0+;  Kp  | log  p | < t < °o  ) . 

Theorem  8 is  proved  in  Section  10. 

The  method  of  proof  of  Theorem  8 uses  the  notions  of  upper  and 
lower  solutions  of  (1.1).  The  necessary  preliminary  material,  which  follows 
the  lines  of  well  known  results  (see  e.  g.  [ 6]  , [ 19] ) is  collected  in 
Appendix  B.  The  inequality  (2.25)  is  established  by  showing  that  the 
solution  of  the  reduced  problem  (1.6)  is  a lower  solution  of  (1.1) 

on  o < t < 00  and  by  showing  that 

w(t,p)  = (t)  + (g(0)  - co(0))exp(-Kt/p)  + y(t)p|log  p| 

is  an  upper  solution  for  suitably  chosen  K and  y (i.e.  that 
t 

-pw'(t)  < f a(t  - s)F(w(t) , w(s))ds  for  0 < t < °o  , where  w(t)  = g(t) 


on  < t £ 0).  Inequality  (2.  25)  then  follows  from  Proposition  2B,  Appendix  B. 


The  question  arises  whether  the  order  0((jl  |log  ^ |)  in  (2.  26)  is 
best  possible.  In  the  linear  case  it  is  not;  for  if  F(y,z)  = y - z,  one 
can  establish  the  inequality 

Kt 

0o(t)  < <t(t,p)  < <tQ(t)  + (g(0)  - 4>0(0))exp(-— ) + y(t)p 

for  0 < t < oo  and  0 < p.  < jjlq , by  the  method  of  proof  of  Theorem  8. 
In  addition,  one  can  compute  a(p)  and  a ^ in  the  linear  case  by  the 
method  of  Laplace  transforms  and  show  that  <*(p)  - is  precisely 
of  the  order  p.  In  the  general  case,  however,  we  have  been  unable  to 
improve  the  estimate  (2.  25). 
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3.  Proof  of  Theorem  1.  The  classical  Picard  successive  approximations 


(or  the  Banach  fixed  point  theorem)  applied  to  the  integrated  form  of  (1.1), 

(1.  2)  show  that  for  each  fixed  p > 0 there  is  a unique  local  solution 
o(t,p)  existing  and  in  C*  on  some  interval  [0,T]  , T > 0.  To  show 
that  this  solution  can  be  continued  (necessarily  uniquely  in  view  of  the 
assumptions)  to  the  interval  [0,°o  )f  it  suffices  to  establish  the 
inequalities  (2.1)  on  any  interval  on  which  the  solution  cf(t,p)  exists. 

For  then  the  solution  satisfies  a priori  upper  and  lower  bounds,  independent 
of  T,  and  hence  can  be  continued  to  the  interval  [0,°°  ) by  a standard 
result  [ 17]  . 

To  establish  (2.1)  on  any  interval  on  which  <f(t,p)  exists  we 
have  from  (1.1),  (1.  2) 

0 

(3.1)  -pcf'(0,p)  = / a(-s)F(g(0),g(s))ds  . 

-00 

The  integral  clearly  exists  since  a«  L*(0,«>)  and  F(g(0),g(s))  is 
bounded  on  (-°o  0]  by  (H  ) H ).  From  (H  ),  a(-s)  > 0 (-°°  < s < 0) 

r Q B 

and  from  (H„),  (H  ),  F(g(0),g(s))  > 0 (-°°  < s < 0),  with  the  strict 

r g 

inequality  holding  for  large  negative  s;  therefore,  4>'(0,p.)  < 0.  Since 
4>  e C1,  one  has  by  continuity  that  <t>'(t,p)  < 0 (0  < t < a),  for  some  a > 0. 
We  claim  first  that 

(3.  2)  4>(t,p)  >1  (0  < t < a)  . 

Indeed,  <t>(0,p)  = g(0)  > 1,  and  by  continuity  (3.  2)  holds  at  least  on 
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some  interval  to  the  right  of  t = 0.  Suppose  0 < t^  < a is  the  first 

point  at  which  <j>(tj,p)  - 1,  and  1 < <^(t , p)  < g(0)  (0  < t < tj) . From  (1.1) 

we  have 

0 ll 

(3.3)  -p<t>  (t.,p)  = / a(t.  - s)F(l,g(s))ds  + / a(t.  - s)F(l,  <t(s, p))ds  . 

1 -oo  1 0 1 

By  (H  ),  a(t.  - s)  > 0 (-»  < s < t ),  and  by  (H  ),  (H  ),  F(l,g(s))  < 0 

a 1 I r Cf 

on  (-00  , o]  , with  the  strict  inequality  holding  near  zero.  Since  <(?  is 
strictly  decreasing  on  [ 0,t^]  , we  also  have  F (1 , 4>(s , p))  < 0 (0  < s < t^). 
Therefore  each  integral  in  (3.3)  is  negative  and  4 1 (t^ , |jl)  > 0,  which, 
in  view  of  4/  (t,  p)  < 0 (0  < t < a) , is  impossible;  this  proves  (3.  2). 

We  next  claim: 

(3.4)  4>'(t,p)  < 0 , 

for  as  long  as  the  solution  exists.  Indeed,  suppose  for  contradiction 
that  a > 0 is  the  first  point  at  which 

(3.5)  <t>'(a,p)  = 0 and  <t>'(t,p)<0  (0  < t < a)  . 

By  the  argument  of  the  preceding  paragraph  we  have 

(3.6)  1 < ^(a, p)  < g(0)  . 

To  prove  (3.4)  we  compute  4>"(ar,p)  from  (1.1)  and  we  shall  obtain  an 
obvious  contradiction  of  (3.5)  by  showing  that 

(3.7)  <j>  "(a,  p)  < 0 ; 

this  implies  that  no  such  a > 0 satisfying  (3.  5)  exists  and  proves  (3.4). 
Indeed,  differentiating  (1.1)  (justified  by  (H  ),  (H  ),  (H  ) - note  that  by 

d 1 Q 
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(H  ),  a'  t L (0,*>  ))  one  has 
a 


r 


(3.8) 


-H-o " (t , ja)  / a‘(t  - s )F (4>(t , hl) , 4>(s,  p))ds 


b 

+ 4>'(t,fi)  / aft  - s)F  (4>(t, (jl) , <})(s,HL))ds  . 


Putting  t = a and  using  (3.5)  gives 

a 

(3.9)  |x)  = / a'(a  - s)F(<J>(or,  jjl) , 4=(s , p))ds  . 

- 00 

Thus  to  prove  (3.7)  we  wish  to  show  that 

a 

(3.10)  1(a)  = J a 1 (a  - s)F(4>(a,  p) , 4>(s , ml) )ds  > 0 . 

_ oo 

We  shall  need  to  consider  two  cases: 

in  case  (i)  a(t)  satisfies  (H  ) with  a(t)  ^ Ae  , A > 0,  X.  > 0 (i.  e. 

a 

a'(t)/a(t)  & -M; 

in  case  (ii)  aft)  satisfies  (H  ) with  a(t)  h Ae  A > 0,  \ > 0. 

a 

Case  (i).  Define  a number  -(3,  (3  > 0,  by  the  relation  4>(a,p)  = g(-(3);  -P 
exists  in  view  of  (3. 6)  and  (Hg).  Since  g may  take  the  constant  value 
4>(ff,n)  on  some  interval  J C (-°°,0],  we  define  -p  uniquely  by  taking 
it  to  be  the  right-hand  end  point  in  such  a case.  We  then  have 
-p  0 

(3.11)  1(a)  = J a'(a  - s)F(g(-p),g(s))ds  + / a' (a  - s)F(g(-p) , g(s))ds 

-oo  - p 

a 

+ f a'fa  - sJFf^fa,  p) , <j>(s , p))ds  . 

0 
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■M 


Since  $'(0,^)  0 we  also  have  from  (1.1) 

-P  0 

(3.12)  0 = J a(o-  - s)F(g(-P),g(s))ds  + f a(a  - s)F(g(-P),  g(s))ds 

_oo  _P 

a 

+ / a(a  - s)F(4>(«,p),<Ms,p))ds  • 

0 

We  next  define  the  function  a by  the  relation 

<r(s)  - a'ja-~  S)  (-«  < s < 0)  , 
a (a  - s)  - ’ 

and  we  observe  that  the  log  convexity  of  a implies  that  <r(s)  is 
negative  and  nonincreasing;  moreover,  since  a'(t)/a(t)  ^ \ > 0,  a(s) 

is  strictly  decreasing,  at  least  on  some  interval  contained  in  (-<»  ,0]  . 
We  rewrite  (3.11)  in  the  equivalent  form 


(3.13) 

where 


- (3  a 

1(a)  = J cr(s)h(s)ds  + J o-(s)h(s)ds 

-P 


-00 


h(s)  = 


a(a  - s)F(g(-p),g(s)) 


(-00  < s < 0) 


^a(a  - s)F(<j>(ar,^),4>(s,|j.))  (0  < s < a)  ; 

we  also  write  (3.12)  in  the  equivalent  form 


(3.14) 


-P  a 

0 = f h(s)ds  + f h(s)ds  . 
-00  - p 


From  the  definition  of  -p  and  (3.6)  one  has  1<  g( - p)  < g (0).  Therefore, 

(H  ),  (K  ) (H  ) imply  that 
a r g 


(3.15) 


h(s)  > 0 


(-00  < s < -p) 
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with  strict  inequality  for  large  negative  s,  and 

(3.16)  h(s)  <0  (-p  < s < a)  . 

Combining  (3.13),  (3.14)  yields 

-(3  a 

(3.17)  1(a)  = J (<t(s)  - n-(-p))h(s)ds  + f (t(s)  - «r( - p) )h(s )ds  . 

-00  _P 

But 

(3.18)  cr(s)  > cr(-p)  (-00  < s < -p)  , 

(3.19)  <r(s)  £ cr(  - p)  (-p  < s < a)  , 

with  strict  inequalities  holding  either  for  s negative  and  large  or  for 
s near  a.  Using  (3.15),  (3.16)  and  (3.18),  (3.19)  in  (3.17)  shows  that 
each  integral  in  (3.17)  is  nonnegative  and  at  least  one  of  them  is  positive. 
This  proves  (3.10)  and  hence  also  (3.7)  and  (3.4).  The  proof  of  the  global 
existence  and  uniqueness  of  the  solution  <j>(t,jx)  of  (1.1),  (1.  2)  and 
of  property  (2.1)  in  case  (i)  is  then  completed  by  a straightforward 
continuation  argument. 

Case  (ii).  The  above  balancing  argument  establishing  (3.  7),  and  hence 
(3.4),  cannot  be  used  when  a(t)  = Ae  >"t,  A > 0,  X > 0,  since  for  this 
case  it  is  readily  established  from  (1.1)  and  (3.9)  that  <(■"(<*,  n)  = 0 
whenever  4j'(ar,|i)  = 0.  Instead  we  proceed  as  follows. 

Define  the  number  y , 1 < yQ  < g(0),  by  the  equation 

/ e F(y0,g(s))ds  = 0.  (For  the  proof  of  the  existence  of  a unique  y 

-oo  U 
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with  this  property  in  the  general  case  of  a(t),  which  also  applies  here, 
see  the  proof  of  (2.  5)  in  Section  4. ) Since  as  in  case  (i)  4'(0,p)<0, 

a repetition  of  the  argument  (3.1)  - (3.  3)  above  with  a(t)  = Ae  , A > 0, 

X > 0,  and  "1"  replaced  by  yQ  in  (3.2)  and  the  paragraph  following 
(3.  2),  shows  that  4(t,|x)  > y for  t>0,  so  long  as  <f'(t,p)<0.  To 

show  that  4 ' (t, m.)  < 0 for  all  t > 0,  differentiation  of  (1.1)  with 
a(t)  ^ Ae‘Kt,  A > 0,  X > 0,  yields  the  equation 

p4'’(t,p)  + 4* (t, |j.)[  Xp  + / e K(t  SV^(<(>(t, p), 4(s, p.) )ds]  = 0 , 

- 00 

1 r ^ Xs 

with  c(0,p)  g(0),  4 ’ (0 , p. ) - - ~ J e F(g(0),g(s))ds  < 0.  Since  > 0, 

^ -00 

X > 0,  p > 0,  an  elementary  argument  shows  that  4 ’ (t , p)  < 0 so  long 
as  the  solution  <fc(t,p)  exists.  This  completes  the  proof  of  (2. 1' ) in 
case  (ii). 

To  prove  (2.2)  subtract  the  equations'  (1.1)  for  4^  and  4.,  and 
apply  the  mean  value  theorem  obtaining 

t 

- P f 4 ’ (t, p)  - 4;(t,p))  = (4,(t,p)  - 4>Jt,p))  J a(t  - s)F  (cr(t), 4,(s, p))ds 

t 

+ J a(t  - s)F2(42(t,p),  t(s))(41(s,p)  - 42(s,p))ds  , 

-00 

where  o-(t)  is  between  4^(t,p)  and  4^(t,p)  and  t(s)  is  between 
4 (s,p)  and  42(s,p).  Let  t = tQ  be  the  last  point  for  which 
4x(t,p)  > 42(t, p);  i.e.  41<t(),p)  = 42(t0,p) , while  ^(t, p)  > 42<t, p) 
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for  -°°  < t < t 


Since  g^tl^g^t),  -oo<t<0,  it  is  clear  that 

t >0.  Both  integrals  exist  since  a c L (0,<»),  F c C , and  4.,  4'? 

0 — 1 c. 

satisfy  (2. 1)  (or  (2. 1')).  From  the  definition  of  t and  F2<0  one  has 
-plcjUg,  p)  - - °>  with  the  ec!uality  sign  holdini3  ^ and  °niy 

if  g1  = g^  on  ( -°o  , 0] . Since  p>0  this  implies  4j(tQ , p)  - 4 > 0 
for  g 3?  gz  and  this  is  impossible. 

To  prove  (2.  3)  put  z(t)  = *(t, Pj)  - 4>(t, P2> • Then  from  (1. 1),  (1.  2) 
one  has  z(0)  0 and 

1 1 r° 

z ' (0)  = (--=-+  — ) J a(-s)F(g(0),  g(sj)ds  > 0 . 

P1  ^2  -00 

By  continuity  suppose  that  there  exists  T > 0 such  that  z(t)  > 0 (0  < t < T) 
and  z(T)  = 0.  Then 

T 

z'(T)  = - — / a(T  - slFWT^^lSjiijNds 

F]  .00 

1 T 

+ — f a(T  - s)F(4(T  p ) , 4(s, p ) )ds  . 

UL  " C C, 

^2 

T 

But  by  (2.1)  / a(T  - s)F(4.(T,p.),4(s,Pi))ds  > 0,  (i  = 1,2).  Therefore 

-OO 

^1  > ^2  fmp^es  ^at 

. T 

z'(T)  >—  / a(T  - s){F(4(T,p2),4(s,p2))  - F(4(T,  pj) , 4(s , p^jds  , 

.00 

which  by  the  definition  of  T and  <}>(*>  M^)  = = g^  on  yields 
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1 r 

(3.20)  z'(T)  > J a(T  - s){F(<j.(T,p  ) , 4 ( s , jx  J)  - F(c(T,  p ) , <J>( s , p ) ) }ds  . 

M-j  Q c*  c.  £ L 

Applying  the  mean  value  theorem  in  (3.  20)  gives 

1 T 

(3.21)  z'(T)  > J a(T  - s)F  U(T,p  ),;(s))(<Ms,p  ) - «|>(s,p.))ds  , 

Pj  0 i.  i.  i.  l 

when  t.(si  is  between  <t(s,p.,)  and  ct ( s , ) . Since  a(T  - s)  > 0 

and  F < 0 by  assumption  and  since  <t(  s,p  ) - <t(s,p  ) <0  on  0<s<T 

^ L.  1 

by  the  definition  of  T,  (3.21)  implies  that  z'(T)  > 0.  Therefore  there 
cannot  exist  a T > 0 such  that  z(T)  = 0;  this  proves  (2.3),  and  completes 
the  proof  of  Theorem  1. 
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4.  Proof  of  Theorem  2 and  Corollary  2.1.  We  observe  first  that  if  ^ 
is  a (continuous)  solution  of  (1.6)  for  t > 0,  then  <fg(0)  = 4>0(°  ) 
must  satisfy 

0 

(4.1)  0 =/  a(-s)F(<f0(0),g(s))ds  . 

- 00 

This  integral  clearly  exists  for  any  number  $ (0)  since  a«  1^(0, °°) 

and  (H  ).  (H  ) hold;  moreover  by  (H  ),  the  integral  in  (4.1)  is  a continuous, 
F ’ g * 

strictly  increasing  function  of  (the  parameter)  bg(0).  The  integral  is 
negative  if  one  chooses  <fcQ(0)  = 1;  it  is  positive  if  one  chooses 
$ (0)  = g(0).  Therefore,  there  exists  a unique  number  <^(0)  = yQ  for 
which  the  integral  in  (4.1)  vanishes  and  clearly  1 < yQ  < 4^(0).  71113 

proves  (2.  5). 

We  next  observe  that  if  <b  is  a continuous  solution  of  (1.6)  on 
0<t<T<°°,T>0,  for  which  <f>  ( 0 ) = yQ,  then  4>Q  is  unique. 

For  if  u,  v e C[  0,T]  , u(0)  = v(0)  = yQ,  u(t)  = v(t)  = g(t)  on  -«><t<0, 
and  u,  v satisfy  (1. 6)  on  [ 0 , T]  , then  by  (H  ) and  the  mean  value 
theorem  we  have 

0 

(4.2)  0 = (u(t)  - v(t))[  / a(t  - s)F  (|  (t),g(s))ds 


I 

J a(t  - s)F1(|2(t),Tll(s))ds] 


t 

+ / a(t  - s)F  (£  (t) , iq  (s))(u(s)  - v(s))ds,  0<t<T, 


for  some  c,  (t),  ^ (t),  ^(t)  between  u(t)  and  v(t),  and  for  some 
r)^[  s ) , ^ (s)  between  u(s)  and  v(s).  But  then  at  L(0,°°),  the 
continuity  of  F,  F , F]  > 0,  and  the  continuity  of  u,  v on  [0,T], 
together  with  Gronwall's  inequality  applied  to  (4.  2),  imply  that 
u(ti  = v(t) , 0 < t < T. 

We  next  deduce  existence  of  solutions  of  (1. 6)  on  [ 0,  T]  , T > 0 

00 

arbitrary,  as  follows.  Let  0<  £ £ t £ T < 00  and  let  {p-n)n  q be 

an  arbitrary  sequence  with  p < p , lim  p^  = 0,  and  consider 

n n-~  00 

GO 

the  sequence  (<t(t,p  ) ) of  solutions  of  the  initial  value  problem 
7 n 1 n 0 

(1.1),  (1.2)  on  e£t£T.  By  Theorem  1,  1 < <|>(t,p  ) < g(0),  and 

c(t,p  ) < c(t,p  ),  e £ t £ T,  n -0,1,  . . . ; therefore,  z(t)  = lim  «j>(t,pn), 
n n n 00 

e<t<  T,  exists.  Moreover,  from 

) z(t  + h)  - z(t)  I £ | z (t  + h)  - 4>(t  + h,pn)  I 

+ | <b(t  + h,p^)  - <t(t,  p^)  I + 1 4>(t , Pn>  * z(t)  I 

for  t,  t + h € [ e,T] , and  from  the  continuity  of  4>(t,p),  it  follows 

that  the  limit  function  z is  continuous  on  [ e,T]  . In  fact,  by  Dini's 

theorem,  the  sequence  {b(t,pn)}  converges  to  z(t)  uniformly  on  [ e,T]  . 

We  also  observe  that  lim  z(e)  exists.  The  next  task  is  to  show  that 

e — 0 

z(t),  with  z(0)  - lim  ^ z(e),  satisfies  equation  (1. 6)  on  0£t£T. 
e — 0 

It  is  easier  to  do  this  by  using  the  estimate  (2.16)  of  Theorem  5 for 
-p<t>'(t,p),  0 £ t £ T;  we  shall  only  indicate  below  how  to  carry  out  the 
passage  to  the  limit  in  (1.1)  without  the  estimate  (2.16).  Consider  the 
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equation  (1.1)  with  p replaced  by  p^,  y replaced  by  c(t,p  ).  We  have 


(4.  3) 


-Pn4>'(t,pn)  J a(t  - s)F(«t>(t, M-n *, «*>(s , jo.  ))ds,  0<t<T, 


cdt,  |i  ) g(t)  (-00  < t < 0) 
n — 


from  (2.16),  lim  p o' (t , p ) 0,  0<t<T.  Ontherighthandsideof 


n — so 


(4.  3)  we  use  at  L (0,*>  ),  the  continuity  of  F,  and  the  uniform 

convergence  of  {b(t  p } to  conclude  that 

n 


lim  J a(t  - s»F(o(t,p  ),b(s,p  )ds  - J a(t  - s )F(z(t) , z(s))ds  (0  < t < T)  . 

nri  * * n 


n — » -ao 


This  shows  that  z(t)  satisfies  the  reduced  equation  (1.6)  on  0 < t < T. 
By  the  uniqueness  of  continuous  solutions  of  (1.6)  we  identify  z(t) 
with  bQ(t).  One  can  avoid  using  the  estimate  (2.16)  by  integrating  (4.3) 
from  0 to  t,  using  the  boundedness  of  the  functions  c(t,p  ) (by 
Theorem  1),  and  by  passing  to  the  limit  in  the  integrated  form  of  the 
equation,  justified  by  the  assumptions;  we  omit  the  details.  Since  T > 0 
is  arbitrary,  this  completes  the  existence  and  uniqueness  of  a continuous 
solution  of  (1.6)  on  0 < t < ao. 

To  establish  the  remaining  conclusions  of  Theorem  2,  we  first 
assume  that  b t C[0,»).  Then  differentiation  of  (1 . 6) , justified  by 
(H  ),  (H  ),  (H  ),  yields 

3 r g 

r*  S 

(4.4)  0 J a'(t  - s)F(b0(t),60(s))ds  + bQ(t)  J a(t  - s)F  (4-  (t),  b (s))ds  , 

— 00  _ qq 
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J 


1 


— 


for  0 < t < go  , where  it  is  useful  to  note  that  the  coefficient  of  c ^ (t ) 
in  (4.4)  is  strictly  positive.  It  is  evident  that  if  cf)  is  a solution  of 
(4. 4)  on  0 < t < T,  satisfying  the  initial  condition 

g(t)  (-<*  < t < 0) 


(4.  5) 


*o(t) 


0 


ft  0)  , 


where  yQ  is  defined  below  (4.1),  then  will  be  a continuous 

solution  of  the  reduced  equation  (1.6)  on  [0,T),  and  by  uniqueness  of 
continuous  solution  of  (1.6)  it  will  be  the  only  such  solution  of  (1.6).  To 
establish  the  complete  equivalence  of  the  reduced  equation  (1.6)  and 
the  initial  value  problem  (4.4),  (4.  5),  one  has  to  show  that  any  continuous 
solution  oQ  of  (1.6)  is  differentiable.  This  can  be  accomplished 
directly  by  forming  difference  quotients  in  (1.6),  applying  the  mean  value 
theorem  and  showing  that  one  can  pass  to  the  limit  thus  arriving  again 
at  (4.4).  Since  for  any  continuous  solution  of  (1.6)  each  integral  in  (4.4) 
defines  a continuous  function  under  our  assumptions,  and  since  the 
coefficient  of  f t ) is  bounded  away  from  zero,  it  follows  that  g^ft) 
is  continuous.  This  completes  the  discussion  of  the  equivalence  of  (1.6) 
and  (4.4),  (4.  5)  which  enables  us  to  deduce  the  remaining  conclusions  of 
Theorem  2 concerning  solutions  of  (1.6)  by  studying  (4.4),  (4.5). 

-\t 

A simple  calculation  shows  that  in  the  special  case  aft)  h Ae  , 

A > 0,  X > 0,  we  have  4>^(t)  = 0 for  t > 0 (from  (4.  4)1,  so  that 
<j>  (t)  = yQ  for  t > 0,  and  this  proves  the  remark  below  (2.6). 
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If  aft)  ^ Ae  , A > 0,  \ > 0,  we  wish  to  show  that  (2.6)  is 
satisfied.  From  (4.4}  we  first  have 

0 

J a'(-s)F(y0,  g(s))ds 
(4.6)  -<J)' (0)  — 

V ' 0 

J a(-s)F  (y  g(s))ds 

-00  1 u 

By  (H^),  (Hf),  (Hg)  the  denominator  in  (4.6)  is  positive.  To  determine 
the  sign  of  the  numerator  we  note  first  that  by  the  definition  of  yQ 

0 

/4‘  n J a(-s)F(y  g(s))ds  = 0 , 

- 00  u 

and  by  the  log  convexity  of  a,  repeating  the  argument  of  Theorem  1 
showing  that  the  integral  1(a)  > 0 in  (3. 10)  (here  we  take  a = 0, 
dKn)  - yQ,  b(s,p.)  = g(s),  and  define  -p,  p > 0,  by  the  relation 
g(-p)  r yQ),  we  find 

0 

^4-81  J a'(-s)F(y  g(s))ds  > 0 . 

- 00  U 

Therefore,  by  (4.6),  (4.7),  (4.8),  one  has  ^(0)  < 0. 

Since  € C (0,<»)  by  continuity  one  has  4>^(t)  <0,  0 < t < a, 

for  some  a > 0.  We  claim,  as  in  the  analogous  part  of  the  proof  of  Theorem  1 
(see  the  proof  of  (3.2)),  that 

f4-9>  bQ(t)  >1  (0  < t < a)  . 

Since  6Q(0)  yQ  > 1 one  has  by  continuity  that  (4.9)  holds  on  some 


-26- 


interval  to  the  right  of  t 0.  Letting  0 < tJ  < a be  the  first  point 
at  which  $ (t  > 1,  and  using  yQ  > cQ(t)  >1  or  0 < t < t^,  one 

computes  c'(t^)  ^rom  <4.4),  and  one  finds  easily  that  ^(t^)  > which 
is  impossible.  This  proves  (4.9). 

We  next  show  that 

(4. 10)  cQ(t)  <0  (0  < t < *>  ) . 

Here  the  procedure  differs  somewhat  from  the  analogous  part  of  the  proof 
of  Theorem  1 in  that  we  do  not  need  to  compute  <?Q.  Suppose  t - a is 
the  first  point  at  which 

(4.11)  o n(a)  0 and  o^t)  < 0 (0  < t < a)  . 


Moreover,  by  (4.4)  we  have 


(4.12) 


0 a 

I a'(a  - s)F(c^(a'),  g(s))ds  + J a (a  - s)F(4>0(ar),4>0(s))ds 


/ a (a  - s)F  (c()  (a),g(s))ds  + J a(a  - s)F  (<t>Q  (or) , <}:  ( s ) )ds 
-oo  0 


The  denominator  in  (4.12)  is  clearly  positive  by  (H  ),  (H  ),  (H  ).  The 

a I Q 

balancing  argument,  involving  (1.6)  in  the  t -a  and  the  log  convexity  of 
a which  has  been  previously  employed  in  the  proof  of  (3.10)  and  (4.8), 
shows  that 


0 a 

J a'(«  - s )F ( <+>0 ( or) , g(s))ds  + J a'(«  - s)F($0(»),<j>0(s))ds  > 0 . 

- oo  0 


Therefore  (b^ta)  < contradicting  (4.11)  and  proving  (4.10).  This 
completes  the  proof  of  (2.6). 
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The  proof  of  property  (2.7)  is  similar  to  that  of  (2.2)  in  Theorem  1 


r 


and  is  omitted. 

To  prove  property  (2.8)  we  observe  first  that  since  <p(0,p)  g(0) 

and  since  <fc0<0 ) yfj,  ( 2.  8 ) is  true  at  t 0 by  ( 2.  5) , and  therefore , 
by  continuity,  (2.8)  is  true  in  some  interval  for  t > 0.  Suppose  that 
o(t,p)  - is  zero  for  the  first  time  at  t tQ  > 0 and  o(t,p)  - <?Q(tl  > 0 

(0<  t < t ).  Then  from  a. 1),  (1.6)  at  t tQ,  one  obtains  by  subtraction 

*0 

(4.13)  -po'(t0,p)  = J a(tQ  - s)[F(<t)(t0,p),4.(s,p))  - F(d0(t0),o0(s))]ds  . 

Applying  the  mean  value  theorem  to  the  difference  under  the  integral 
(note  that  <t>(t0 , h->  = ^0(t0)  by  the  definition  of  tQ)  and  using  F2<y,z)<0 

and  ®(s,p)  - <ts  ( s ) > 0 (0  < s < tQ),  it  is  evident  that  the  right-hand 
side  of  (4. 13)  is  negative.  This  implies  that  d'(t0,p)>0  which 
contradicts  (2.1)  and  proves  (2.8);  this  completes  the  proof  of  Theorem  2. 

To  establish  the  conclusion  (2.10)  of  Corollary  2.1  let  n > 0 be 
given.  Using  (2.4),  (2.8),  (2.9)  choose  T > 0 so  large  that 

|a(p)  - <f(t,u)  I < ^ and  0 < <•>  (t)  - * < for  t > T.  By  the  convergence 

of  d>(t,u)  to  $ (t)  as  p -*  o’1  on  0 < t < » choose  p>0  sufficiently 

small  that  0 < <f(T,p)  - 4>Q(T)  < Then 

o < a(p)  - < |a(p)  - 4>(T,  p)  I + (b(T,p)  - oQ(T))  + UQ(T)  - <*Q)  < n , 

for  T > 0 sufficiently  large  and  p > 0 sufficiently  small,  proving  (2.10). 
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(5.1) 


g(tl 


5.  Proof  of  Theorem  j.  In  view  of  the  monotonicity  property  (2.7)  of 
solutions  of  (1.6)  with  respect  to  the  function  g,  it  suffices  to  prove 
the  result  for  the  function  g given  by 

1+6  if  -g  < t < 0 
1 if  t < -t],  f > 0,  n > 0 . 

Note  that  any  such  discontinuous  function  g can  be  approximated  by 
arbitrarily  smooth  functions  g satisfying  assumption  (H  ).  Since 
1 < g(t)  < 1 + 6,  -oo  < t < 0,  property  (2. 6)  implies  that  1 < <t>  (t)  £1+6, 
0 < t < oo  . This  means  that  the  arguments  of  F in  (1.6)  are  close  to  1, 
if  6 > 0 is  sufficiently  small  which  will  be  the  case  in  what  follows. 
For  this  reason  we  assume  consistent  with  (H  ) and  without  loss  of 
generality  that 

(5.  2)  F1(l,l)  1,  F2(l,l)  = -1  ; 


note  that  F(x,x)  = 0 (x  > 0)  implies  that  Fj(x,x)  = -F  (x,x). 
Substituting  (5.1)  into  (1.6)  yields  (note  a t L (0,oo), 


-r\  y 

(5.3)  F(<t  (t),l)  / a(t  - s)ds  + F(<ts  (t) , 1 + 6)  / aft  - s)ds 


- 00 


J aft  - s)F(*  (t),4>  (s))ds  = 0 . 


Using  (5.  2),  the  mean  value  theorem,  and  F € C (IR  x ]R)  yields 
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(5.4) 


(o  (t)  - 1 J a(t  - s)ds  + (4>0 (t)  - 1 - 6)  J alt  - s)ds 


0 t 

o(6)  j a(t  - s )ds  + / a(t  - s)($  (t)  - 4>Q(s))ds 
- jo  0 


+ / a(t  - s)[o((p0(t)  - <j>0(s))]ds  0 , 

0 

where  (by  a « L1^,'*’),  aft)  > 0)  the  terminology  w(t)  = o(6)  / a(t  - s)ds 

-00 

00 

means  that  for  every  e > 0 one  has  |w(t'  I < e6  J a(4)d|  for  t > 0 

t 

and  for  6>0  sufficiently  small.  An  equivalent  form  of  (5.  4)  is 

00  t t+T) 

(5.5)  (d  (t)  - 1)  f a(fe)d|  - / a(t  - s)(<(>0(s)  - l)ds  = 6 / a(i)d| 


+ o(6)  / a(£)d|  + / aft  - s)[ o(«t>  (t)  - «t>Q(s))]ds  . 
t 0 

Since  1 < <;0(t)  < 1 + 6 on  0<t<oo(  (Hp)  implies  that  the 
first  term  in  (5.  3)  is  positive,  while  the  second  and  third  terms  are 
negative.  Thus  the  third  term  in  (5.  3)  must  be  in  modulus  less  than  the 

first  term.  This  in  turn  implies  that  the  last  term  in  (5.  5)  (or  the  last 

00 

integral  in  (5.4))  must  be  o(6)  / a(|)d|.  Therefore,  putting  z(t)  = <t>Q(t)  - 1 

t 

in  (5.  5)  yields  the  equivalent  equation 


(5.  6)  z(t)A  - a * z(t)  = 6 / a(s)ds  + o(6  f a(s)ds)  (0  < t < °°)  , 

t t 

00 

where  * denotes  the  convolution  and  A = J a(s)ds. 

0 


■■■■■■■■I 


We  shall  now  apply  a simple  Tauberian  theorem  for  (real)  Laplace 
transforms  [ 20;  Theorem  4.  3,  p.  192]  to  solutions  of  (5.6).  Put 


t-h 


(5.  7 i 


iMt)  = J a(s)ds,  co(t)  = f a(s)ds  , 
t t 


(5.8) 


z (p ) = / e Ptz(t)dt 
0 


(p  > 0)  . 


.1. 


Since  z is  bounded  and  continuous  on  [ 0 , 00 ) , and  since  a « L (0,°°),  z(p) 
as  well  as  o(p)  and  w(p)  exist  for  p > 0.  Noting  that  A = a(0)  and 
that  multiplication  of  (5.  6)  by  e~pt  and  integration  preserves  the  o 

A 

relation  when  p is  real,  we  obtain  on  solving  for  z(p) 

(5.9) 


;«>»  ■ T~r ,p  > oi 


By  assumptions  (H  , (2.11),  and  Lebesgue's  dominated  convergence 
theorem,  and  integration  by  parts  we  also  have 

(5.10) 


w 


(0)  = f J a(s)dsdt  = J sa(s)ds  , 

0 t 0 

^ 00  t+T]  00  00  00 

(5.11)  ^j(0)  = f f a(s)dsdt  < f J a(s)dsdt  = f sa(s)ds 

0 t 0 t 0 

Moreover,  by  Fubini's  theorem  and  (2.11) 


(5.12) 


~~  a(p)  = - J te  Pta(t)dt  (p  > 0)  , 
aP  o 


so  that  by  the  mean  value  theorem  and  (2.11) 
(5.13) 


00 


a(0)  - a(p)  ~ p J ta(t)dt  (p  — 0 
0 
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1 


Hence  using  (5. 10),  (5. 11),  (5. 13),  letting  p - o'*  in  (5. 9),  and  taking 
6 > 0 sufficiently  small,  there  exists  a constant  K(6)  > 0 such  that 

(5.14)  z(p)  ~ ^ (p  - 0+)  . 

The  above  mentioned  Tauberian  theorem  then  implies 

(5.15)  z(t)  ~ K(6)  > 0 (t  - + »)  , 
and  by  the  definition  of  z (5.15)  yields 

<fcQ(t)  ~ 1 + K(6  ) > 1 (t  -*  + oo)  , 

for  6 > 0 sufficiently  small.  This  completes  the  proof  of  Theorem  3. 
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6.  Proof  of  Theorem  4.  Using  (2.13),  equation  (1.6)  may  be  written  in 
the  form 


0 


(6.1) 


Ay (t)  - a * y(t)  J a(t  - s ig(s)ds  (0  < t < f-»)  , 


CO  t 

where  A J a(s)ds  and  a * y(t)  J aft  - s)y(s)ds.  Letting,  as  in 
0 0 

the  proof  of  Theorem  3,  y(t)  z(t)  - 1, 

(6.1)  becomes  the  linear  Volterra  equation 

(6. 2)  Az(t)  - a * z(t)  G(t)  (0  < t < n)  , 

where 


0 


(6.3) 


G(t)  j a(t  - s)(g(s)  - lids 


(0  < t < --a)  . 


We  proceed  as  in  the  proof  of  Theorem  3,  but  taking  any  g t (H  ), 

0 

rather  than  the  special  g of  (5.1).  Taking  the  Laplace  transform  and 
noting  that  A a(0)  we  obtain 

-giej 


(6.4) 


z(p)  = — 


a(0)  - a (p ) 


By  Theorem  2 lim  z(t)  z exists  and  z >0.  If  z >0,  a standard 

GO  30  — 30 

t -*•  00 

Abelian  theorem  for  Laplace  transforms  ( 20,  Cor.  lb,  p.  182]  states  that 
z 

A 00  -f 

z(p)  ~ (p  — 0 ).  We  shall  apply  this  result  to  (6.4)  for  p real. 

Since  by  (H  ) g(t)  > 1 (-°°  < t < 0),  and  since  g is  nondecreasing 
g - - ’ 

with  g(t)  > 1 near  t 0,  we  have  by  Fubini's  theorem 


-33- 


(6.  5) 


uO  0 

0 < G(p)  J e Pt  ) a(t  - s)(g(s)  - l)dsdt 
0 -oo 

0 °o 

J (g( s)  - 1)  J e Pta(t  - s)dtds  (p  > 0) 
- oo  0 

Let  e > 0 be  given;  choose  a number  N N(e)  > 0,  and  using  (H^) 
divide  the  interval  (-°°,0]  into  two  parts,  such  that 

(6.6)  g(s)  - 1 < e (s  < -N)  . 

Then  (6.6)  used  in  (6.  5)  yields 

A 0 00  o 

0 < G(p)  < e j e PS  J e P°a(0)dOds  + 0(1)  (p  — 0 ) , 
-oo  -s 

or  equivalently 

(6.7) 

By  an  integration  by  parts  and  a t 1/(0, 00 ) we  also  have 

00  00 

(6.8)  a(0)  - a(p)  = J (1  - e pt)a(t)dt  = J (1  - e Pt)ePte  pta(t)dt 


00  00 

0 < G(p)  < £ J eP  ' J e ^"a(0)dOd(r  + 0(1)  (p  -*  0 

0 <r 


per  r -p0_  . -4 


.1, 


0 0 

00  oo 

- P J ePt  / e p9a(0)d0dt  >0  (p  > 0)  . 

0 t 

Therefore,  (6.7),  (6.8)  substituted  into  (6.4)  yields 

oo  oo 


(6.9) 


e I ePt  I e P°a(0)dGdt  4 0(1) 
z (p)  < — (p  - 0+)  . 


oo  oo 


p / ePt  J e p6a(0)d0dt 

0 t 
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But  by  (2. 14) 

00  00 

lim  / ePt  / e"P°a(G)d0dt  + * , 
p-0+  0 t 

and  this  fact  used  in  (6.9)  shows  that  for  any  e > 0 

z(p)  < ~ + o(— ) (p  -*■  0 ) , 

- p p 

since  e > 0 is  arbitrary^  ^ 6 is  impossible.  This  completes  the 

proof  of  Theorem  4. 
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1 • Proof  of  Theorem  5 and  Corollary  5.1.  Let  c be  the  solution  of  (1. 1) , 
(1.2).  We  return  to  equation  (3.8)  obtained  by  differentiating  (1.1;  and 


we  write  (3.8)  in  the  form 

(7.1)  - m-«|>  " (t , (x)  G(t,|j0<t,(t,|i)  + f(t,(j.)(0<t<‘»), 

where 


(7.  2) 


G(t,fx)  J a (t  - s)F^(d(t,  |x),  <f(s , m- ) )ds  , 


(7.  3) 


f(t,|x)  = / a 1 (t  - s )F(d>(t, fi.) , cj)( s , (j.) )ds  . 


Since  Fj  > 0,  at  L (0,°°),  a(t)  >0  (0  < t < °°),  and  since  <t  satisfies 


(2.  2),  (2.4),  (2.6),  we  have  by  (H  ), 


(7.4) 
where 

(7.5)  { 


0 < y A < G(t, |x>  < I'A  (0  < t < °°,  fji  > 0)  , 


r y inf  F (y , z ) , F = sup  F(y,z),  S = [ 1,  g(0)]  X [ 1 , g(0 )]  , 
1 S S 


A = J a(s)ds 
^ 0 


We  next  show  that  there  exists  a constant  K > 0,  independent  of  p, 
such  that 


(7.6) 


I f (t , )jl ) | < K f a(s)ds  (0  < t < fx  > 0 ) . 

t 


Write  f(t,p)  = Ij  + I , where 
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0 

Ij  j a'(t  - s)F(<*>(t , p) , g(s ))ds 

- JO 

t 

I a * ( t - s)Ffo(t,  (j.) , ( s , p))ds  . 

0 

Consider  I,  first  and  recall,  from  the  proof  of  Theorem  1,  that 

F(<t(t,  (jl)  , <t  (s , p) ) <0  (0  < s < t).  Since  a'  < 0,  one  has  on  letting 

>(t  - si  ,Q  ~ S~  and  on  using  the  log  convexity  of  a,  that 
alt  - s)  ’ 

t 

(7.7)  0 • I < ( - " (0)  i ! J a(t-s)F(c{:(t,p),4(s,p))ds|. 

L 0 

But  now  (1.1),  conclusion  (2.1)  of  Theorem  1,  together  with  a simple 
consideration  of  signs  of  the  terms  on  the  right-hand  side  of  (1.1)  shows  that 
t 0 

(7.8)  | J a(t  - s)F(o(t, p) , o(s , p))ds  | < J a(t  - s)F(c (t , p) , g(s) )ds  . 

0 -oo 

Moreover,  the  boundedness  of  c,  g and  the  continuity  of  F imply  that 

0 oo 

(7.9)  / a (t  - s)F(t|  (t,p.),g(s))ds  < sup  I F(e  (t,  p),  g(s))  | J a(s)ds  . 

-33  0<t<°c  t 

-oo  <s  < 0 

Combining  (7.  7) . (7.8).  (7.9)  shows  the  existence  of  a constant  K, 
independent  of  p,  such  that  satisfies  the  estimate  (7.6).  In  a 

similar  way  one  finds  on  using  the  log  convexity  of  a that 

00 

ll  | < (--t(0))  sup  |F(<f;(t,p),g(s))  I / a(s)ds  . 

0 < t < oc  t 

-oo  < S < 0 


Combining  the  estimates  for  1^  and  I establishes  (7.6). 


r 


Returning  to  (7.1)  we  have  by  integration 


(7.10)  ( - o ' (t , jjl ) exp(—  J G(s,p)ds))  exp(~  J G(s,p)dsi. 


By  (7.6)  one  has,  for  0 < t < 00  and  p > 0, 


(7.11) 


K 1 

exp (— J G(s,  p)ds)  < — exp(~  J G(s,  pjds  + log  J a(sids)  . 


P P 


Before  integrating  (7.10)  observe  that  by  a,  a'  e L (O,00)  and  by  the  log 


convexity  of  a one  has 


H 1 r 

— (-  J G(s,p)ds  + log  J a(s)ds) 


G(t.p)  am 


J a(s)ds 


I a (s)ds 

G(G  p)  + 1 ^(t'P.1 

P 00  P 

J a(s)ds 
t 


0(1)  (0  < t < « ^ > 0)  . 


A simple  but  tedious  calculation  then  shows  that  using  (7.4),  (7.6)  in  (7.11) 
there  exist  constants  p^  > 0,  K - K( p^ ) > 0 such  that 

(7.12)  |/  exp(—  f G(s , p)ds)d£,  | < 


1 t « 

K exp(_  J G(s,p)ds  + log(  J a(s)ds))  (0<t<°°;0<p<  p ) . 


Then  integrating  (7.10)  and  using  (7.12),  and 


1 r 

-<t'(0,n)=-J  a(-s)F(g(0),g(s))ds  , 
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as  well  as  conclusion  (2.1)  of  Theorem  1,  yields  the  estimate  (2.16), 


where  \A.  This  completes  the  proof  of  Theorem  5. 

Proof  of  Corollary  5.1.  Integrating  (2.16)  from  t to  infinity  and  using 
conclusions  (2.1)  and  (2.4)  of  Theorem  1 yields 

-K  t/p  ~ °°  00 

0 < 4>(t,fi)  - or((j. ) < Ke  + K J (J  a(s)ds)ds  (0  < p < PQ : 0 < t < 00 ) ; 

t t 

an  integration  by  parts  gives 

~ -KiW  ~ “ 

(7.13)  0 < o(t,p)  - «(p)  < Ke  + K J (s  - t)a(s)ds  (0  < p < pQ;  0 < t < «)  . 

t 

Since  a(t)  is  log  convex,  log  a(t)  is  bounded  below  by  an  affine 

-pt 

function;  thus  there  exist  constants  a > 0,  p > 0 such  that  aft)  > ae 
(0  < t < *>).  Hence  given  any  tQ  > 0 there  exists  a constant  which  we 
again  call  pQ  > 0 such  that  the  integral  term  in  (7.13)  dominates  the 
first  term  for  t < t < °°,  0 < p < pQ ; this  statement  remains  true  at  tQ  - 0 
by  possibly  increasing  K.  This,  together  with  (7.13)  establishes  (2.17). 

The  first  inequality  in  (2.18)  follows  from  (2.8).  To  prove  the  second 
inequality  in  (2.18)  note  that 

g(t,p)  - <i>  (t)  = 4>(t,p)  - a(p)  + o(p)  - o0(t) 

< <)>(t,p)  - «(P)  + «(P)  - «Q  , 

where  the  last  step  follows  by  Theorem  2;  (2.18)  now  follows  by  estimating 
d(t,p)  - «(p)  by  (2.17).  This  completes  the  proof  of  Corollary  5.1. 


-39- 


>mnBOBHnVr  ■•■  Wmmm  — .. 

(8.  Proof  of  Theorem  6.  Let  c.  be  the  unique  solution  of  (1.6)  (see 
Theorem  Z).  Thus  o„  satisfies 

I 

0 t 

(8.1)  J aft  - s)F(<t>  (t),g(s))ds  + J a(t  - s)F(<p  (t) , (j>0(s))ds  0(0<t<^>. 

By  ae  L^(0 , 00 ) , (H^),  (H  ),  and  the  boundedness  of  cQ  there  exists 

a constant  K > 0 such  that 

0 *> 

(8 . Z)  0 < J a(t  - s)F(o0(tl,g(s))ds  < K J a(s)ds  (0  < t < »)  . 

-00  t 

(The  first  inequality  in  (8.  Z)  follows  by  simple  consideration  of  signs  in  (8.1)). 
Applying  (H^)  and  the  mean  value  theorem  we  have 

F(cQ(t),ct0(s))  = F1(<b0(s),o0(s))(c0(t)  - o0(sl)  + o(60(t)  - <Q(s))  , 

which  by  the  boundedness  of  Cq,  F^  > 0,  and  the  continuity  of  F^ 
implies  that  there  exist  constants  > 0 such  that 

(8.3)  M^fs)  - d0(t))<  -F(60(t),<t0(s))  < M2(o0(s)  - d>0(t)):  0 < s < t < =«  . 

Using  (8.Z),  (8.3)  in  (8.1)  shows  that 
t 

(8.4)  f alt  - sild^ft)  - *Q<s))ds  - 4>(t)  (0  < t < 30 ) 

0 


where 

(8.5) 


00 

4>(t)  = Olf  a(s)ds)  . 

t 


k 


-40- 


Let 


(8.6) 


z(t)  c (t)  - «0  («0  lim  oQ(t))  • 

t “*■  00 


Then,  from  (8.4),  z satisfies  the  linear  Voltcrra  equation 


(8.7) 


t t 

z(t)  J a(s)ds  - J aft  - s)z(s)ds  o(t)  (0  < t < *>) 
0 0 


writing  / a(s)ds  / a(s)ds-  / a(s)ds  and  noting  that  z(t)  is 
0 0 t 

bounded  shows  that  (8.7)  may  be  written  in  the  form 


(8.8) 


z(t)A  - a * z(t)  = w(t)  (0  < t < °°;  A f a(s)ds  > 0) 

0 

OC  00 

w(ti  = ^(t)  + z(t)  j a(s)ds  CM  J a(s)ds),  0 < t < 00  . 


We  now  solve  (8.8)  by  Laplace  transforms.  By  the  argument  of 

oo 

Theorem  3 and  using  J ta(t)dt  < 00 , we  find  that  if 

0 

00 

w(0)  J w(t)dt  # 0 
0 

(note  that  this  integral  exists  in  view  of  the  estimate  satisfied  by  w 

and  ta(t)  « L1(0 , °°)),  then  lim  z(t)=z<x>*0.  But,  by  Theorem  2, 

t-  00 

z = 0 . Therefore, 

00 


j w(t)dt  = 0 . 
0 
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M 


Now  integrate  (8.8)  over  [0,T],T>0,  obtaining 

T T T 

A J z(s)ds  - J a * zftidt  J w(t)dt 
0 0 0 

But,  by  (8.9/,  the  last  equation  can  be  written  as 


(8.10) 


T T oo 

A J z(s)ds  - J a * z(t)dt  - / w(t)dt 
0 0 T 


Interchanging  the  order  of  integration  in  the  double  integral  on  the  left- 

hand  side  of  (8.10)  and  using  the  estimate  for  w in  (8.8),  as  well  as 

ta(t)  € l/(0,*>)  yields 

.T  T T oo 

(8.11)  A J z(s)ds  - J z(s)  J a(i)d<rds  0(J  (t  - T)a(t)dt)  (0  < T < <*>)  . 
0 0 0 T 

Using  the  definition  of  A and  combining  the  two  integrals  on  the  left-hand 


side  of  (8 . 11)  yields 


(8.12) 


J z(s )f  a(»)d(rds  0(J'  (t  - T)a(t)dt)  . 
0 T-s  T 


finally,  using  the  fact  that  z is  decreasing  on  0 < T < <*  we  obtai 


from  (8 . 12) 


T T 


J f a('T)d<rds  = 0(f  (t  - T)a(t)dt)  , 
0 T-s  T 

which  on  interchanging  the  order  of  integration  yields 


(8.13) 


z(T)  f o-aMdtr  r Olf  (t  - T)a(t)dt 
0 T 


This  implies  (2.19)  and  completes  the  proof  of  Theorem  6. 
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9_.  Proof  of  Theorem  7.  By  standard  results  the  initial  value  problem  (2.  22) 
has  a unique  local  solution  w £(t)  having 

0 

(9.1)  -4 ' (0 ) J a(-s)F(g(0),g(s))ds  > 0 . 

- 00 

Thus  % decreases  initially.  To  continue  the  local  solution  we  proceed 
similar  to  the  proof  of  Theorem  1.  First,  we  show  that  for  as  long  as 
4 '(t)  < 0 one  has 

(9.2)  &(t)  > yQ  , 

where  y^  is  (see  Theorem  2)  the  unique  value  for  which 

0 

(9.3)  / a(-s)F(y  g(s))ds  = 0 . 

- 00  U 

We  observe  that  equation  (2.  22)  may  be  regarded  as  an  autonomous 
ordinary  differential  equation  having  the  point  y^  as  its  only  critical  point. 
Recall  from  Theorem  2 that  (4(0)  =)g(0)  > yQ.  If  § assumes  the  value 
yQ  at  t Tj,  then  by  (9.2) 

0 

-£'(t1)  = / a(-s)F(y0,g(s))ds  = 0 , 

- oo 

which  is  impossible  if  |'(t^)  < 0.  This  proves  (9.2).  On  the  other  hand, 

4 * (0 ) < 0 implies  by  continuity  that  £'(t)  < 0 for  0 < t < a,  a > 0.  If 
a > 0 is  the  first  point  at  which  4'(a)  = 0,  then 

0 

0 = -^(a)  = J a(-s)F(£(or),g(s))ds  • 

-00 
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This  by  uniqueness  of  implies  that  £(a)  y and  therefore, 
w |(t)  and  w = Y0  would  be  two  solutions  of  (2.  22)  through  the 
point  (“,Y0>,  contradicting  uniqueness.  Therefore, 

(9.4)  4'(t)  < 0 , 

for  as  long  as  the  solution  exists.  Now  (9.2),  (9.4)  and  a standard 
continuation  argument  yield  the  global  existence  and  uniqueness  of  the 
solution  w £(t)  °f  (2.22)  such  that 

(9-5)  4'(t)<0,  yQ  < |(t)  < g(0)  (0  < t < =o)  , 

which  implies  that  lim  £(t)  - L exists  and  £ >y  . 

oo  = oo  — ' o 

T — 00  U 

To  prove  (2.  23)  we  combine  (2.  22)  and  (9. 3)  obtaining 
dt  ° 

(9-6)  ' dr  j a(-s)(F(|  (t),  g(s))  - F(yQ,  g(s))ds  (0  < t < oo)  . 

Applying  the  mean  value  theorem  and  (H^)  yields  the  existence  of  a 0(t), 
0 < 0(t)  < 1,  such  that 

d£  r° 

(9-  7>  - dr  = J a(-s)F1(yQ  + 0(t)(4(t)  - yQ),  g(s))ds[  4(t)  - yQ]  . 

00 

Let  s r f yQ,g(0)]  X [ l,g(0)]  , A / a(s)ds  > O.  By  (H  ) 

0 F 

Y inf  F , (y,  z)  > 0 . 

(y,z)e  s 

Therefore  using  this  and  (9.  5)  in  (9.  7)  gives  the  differential  inequality 
(9.8)  - ^ > yA(£(t)  - yQ)  (0  < t < oo)  . 
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Integrating  (9.8)  and  using  the  initial  condition  £,(0)  g(0)  yields 

the  second  statement  in  (2.  23),  which , together  with  the  existence  of 
lim  |(t)  % also  implies  the  first  statement  in  (2.  23). 

T — * JO 

Applying  the  transformation  (2.  21)  it  is  clear  that  y(t)  §(t/p)  is 
the  unique  solution  of  the  initial  value  problem  (2.  20)  for  p > 0,  0 < t < °c. 
Thus 

r 0 

-p  TT  i(~)  J a( - s )F(|  (— ) , g(s  ))ds  (p  > 0,  0 < t < <») 

(9.9)  < ^ -»  * 

v-  4(0)  = g(0)  . 

Let  <jj(t,p)  be  the  unique  solution  of  (1.1),  (1.  2)  on  0 < t < °°,  p > 0 
(see  Theorem  1),  which  can  be  written  in  the  form 

d r° 

(9.10)  -p  — <|>(t,p)=  J a(-s)F((jj(t  ,p),g(s))ds 

- co 

0 

+ t J a'(-s  + 0 (t , s)t)F(<t(t  ,p),g(s))ds 

- co 

t 

+ f a(t  - s)F(c(t  , p) , <t>(s , p))ds  (0  < t < °o,  p >0)  , 

0 

where  0 < 9(t,s)  < 1 comes  from  the  application  of  the  mean  value 
theorem  to  a(t  - s)  - a(-s).  Not  e that  since  F(4>(t  , p),  g( s ) ) is  bounded 
and  a'  € L^(0,°°),  the  second  integral  in  (9.10)  clearly  exists.  Subtracting 
(9.9)  from  (9.10),  using  the  mean  value  theorem  once  more,  and 

-f 

observing  that  the  last  two  terms  in  (9.10)  are  Oft ) as  t — 0 , 
uniformly  in  p > 0,  yields 
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w 


(9.11) 


-p  - uh]  91(t,p)[  <Mt,p)  - 1(^1]  + O(t) 

(t  • 0 , uniformly  in  p > 0 1 , 

where 

0 

(9.12)  SI  ( t , |x ) J a(-S)F.(|(-)  + 0(t,p)«:(t  ,p)  - £(-)), g(s))ds  , 

Yh  M* 


0 < 0(t, |x)  < 1.  Letting  S [ 1,  g(0l]  x [ 1,  g(0)]  . A j a(s)ds,  and 

0 

using  (H  ),  (H„),  (H  ) implies 
a ’ T g 


(9.13) 


Inf  21  (t,  p'l  > yA  > 0 , 

°<t<t0 
u > 0 


where  y infF.(y,z)  is  independent  of  p.  Using  (9 . 13)  in  (9 .11), 

S 

integrating  (9.11),  and  applying  the  initial  condition  <t>(0,p)  t,(Oi  g(0), 
yields  the  existence  of  a constant  K K(tQ),  independent  of  p such  that 


t - ~L~~  ( t — s ) 
l<t>(t,p)  - £(“)  I < K J ^ e 1 


ds  (0  < t < tfl ; p > 0)  . 


(9.14) 


^l< 

VI 

[*♦*• - J 

(0<  t<  tQ,  p>  0)  , 


£ 

which  is  (2.  24)  with  K = — . This  completes  the  proof  of  Theorem  7. 

yA 
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10.  Proof  of  Theorem  8.  We  show  that 


(10.1)  w(t,p)  - 4>.(t)  + (g(0)  - <t--n(0))exp(-  — ) + |log  ii  | , 

u u fi. 

where  <{)0  is  the  solution  of  (1.6),  is  an  upper  solution  of  (1.1)  for 
0 < t < 30  for  suitable  choices  of  the  constant  K > 0 and  the  function  y . 

By  Appendix  (B)  it  suffices  to  show  that  w defined  by  (10 . 1)  for  0 < t < 30 
and  w(t)  g(t)  (-00  < t < 0)  satisfies  the  integrodifferential  inequality 

t 

(10.2)  -pw'(t,p)  < / a(t  - s)F(w(t,p),w(s,p))ds  (0  < t < 00 ) . 

- 00 

We  begin  by  defining  y and  K in  (10.2).  For  reasons  which 
will  become  apparent  below  let 

t 30 

(10.3)  \(t)  = exp(K  / (/  a(-r)dTidrr)  (0  < t < 30 ) 

0 o- 

where  > 0,  > 0 are  constants  specified  below.  The  function  w 

defined  by  (10.1)  satisfies  the  inequality 

00  00 

(10. 4)  l<w(t,|i)  < d (0)+ g(0)  + exp(K1  J J aMdrdi  )p  | log  p I (0  < t < *)  . 

0 a 

Let  J denote  the  closed  interval  [ 1,  2(4  Q(0)  + g(0))]  and  let  R denote 
the  rectangle  Jx  J.  Let  M > 0 be  a constant  such  that 

(10.5)  max  { 1 F(y , z ) | , I F^(y , z ) | , I F ^(y , z ) | , |F^(y,z)|,  |f  ^(y , z ) | , |l  22<y , z 1 ! < 
(y,z) « R 


and  let  m > 0 be  a constant  such  that  (see  (H  )) 

(10.6)  min  F ,(y,z)  > m . 

(y , z) e R 


Then  for  any  choice  of  \ K independent  of  p,  one  has  for  p 
sufficiently  small  that  the  points: 

(w(t,p),w(s,p)),  (w(t,p),c>0(s)),  (<t>0  (t) , w(  S , (x) ) , (cQ(t)  , Oq(s))  £ R 

for  0 < s < t < : this  statement  is  also  true  for  t < 0 since 

w(t,p)  cQ(t)  g(t>  f°r  t < 0.  Therefore,  the  values  of  F and  its 
first  and  second  partial  derivatives  at  these  points  are  by  (10.  5)  bounded 
by  M . 

00 

1 r 

We  define  K in  (10.1)  by  K ~ rnA,  where  A J a(s)ds. 

0 

With  these  definitions  of  K and  y it  remains  to  verify  (10.  2). 
We  begin  by  doing  this  on  the  interval  0 < t < — p |log  p |.  To  simplify 
the  exposition  let  RHS  denote  the  integral  on  the  right  hand  side  of  (10.2) 
for  t on  any  interval  under  consideration.  We  shall  also  surpress  the 
parameter  p in  w(t,p)  when  no  confusion  arises.  By  the  mean  value 
theorem  we  have  for  0 < t < ^ p |log  p I : 
t 

(10.7)  RHS  J a(t  - s)F((w  - $ )(t)  + <f0(t),(w  - <f>0)(s)  + ci>0(s))ds 

- 00 

t 

= (w  - <fj(t)  f a(t  - s)F,te(t),4(s))ds  + 0(p  |log  p |) 

0 J 1 

- oo 

where  the  O term  does  not  exceed  ^ |i0g  ^ | ; and  £(t)  lies 

between  «j>  (t)  and  w(t).  Thus  (10.1),  (10 . 6)  used  in  (10 . 7)  yields 

(10.8)  RHS  > mA[  (g(0)  - 4>0  ( 0 ) )e " Kt//hX  + \(t)p  |log  p I]  + 0(p  |log  pi)  . 
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Choosing  ^ in  (10.3)  so  that  mAyQ  > 4Ma(0 )/K  and  leaving 
arbitrary  at  the  moment,  (10.8)  becomes 

(10.9)  RHS  > mA[  (g(0 ) - o0(0))e  Kt//,'t  + -^  yQ|i  llog  p.  |]  . 

On  the  other  hand,  it  follows  from  (10.1)  and  (10.  3)  using  y'(t)  > 0 (0  < t < °°) 
that 

(10.10)  -pw'(t)  < -HL^^(t)  + K(g(0)  - c0(0))e'Kt/M"  (t  > 0)  . 

Since  -Cg(t)  > 0 it  follows  from  the  choice  of  K by  comparing  (10.9), 

(10.10)  that  for  p sufficiently  small  the  desired  inequality  (10.  2)  holds 

on  the  interval  0 < t < ^ p |log  p I . 

We  next  verify  (10.  2)  on  the  interval  ~ p [log  p | < t < K where  K 

I\  Z Z 

is  a positive  constant,  independent  of  p,  to  be  determined.  Observe 

that  for  t > — p | log  p | , 

K 

(10. 11)  (w  - <j>  )(t)  1 y(t)p  I log  p I + 0(p  ) . 

Thus  by  the  mean  value  theorem  and  the  fact  that  w(t)  = o^ft)  - g(t)  (-00  < t < 0) 
one  has 

t 

(10.12)  RHS  = (w  - <)>  )(t)  J a(t  - s)F  (6(t),fc(s))ds 

- oo 

2p  | log  p |/K 

+ / aft  - s)(w  - c{:0)(s)F2(^(t),|(s))ds 
0 

t 2 

+ f aft  - s)(w  - <1  )(s)F  ,(|(t),£(s))ds  (t  > - p I log  pi), 

2p  ] log  p l/K  ° 2 K 
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where  t, ( t ) is  between  o()<t)  and  w(t). 

first,  second,  and  third  integral  in  (10.12). 

■,  nondecreasing  one  has 

2Ma(0)  . 2],  1 2 

(10.13)  |I2I<  — £ — l log  I1  I 


Let  I , 1^,  I3  denote  the 
Then  by  (10.1),  (10.5),  and 


+ 


Ma(0 1 


fi(g(0)  - ( o > ) , 


ana  similarly 


'P  S' 

(10.14)  |l  | < M[  y(t)n  | log  M.  I + 0(|JL  )]  J a(t  - s)ds  . 

f llog  m-I 


Note  that  I , I are  each  negative  while  1^  is  positive  and  for  p.  > 0 
sufficiently  small 

(10.15)  I)  > ~ mA{y(t)p.  ilog  ^ ! + 0(p.  )}  . 

So  if  is  chosen  (independent  of  p ) and  sufficiently  small,  it 

follows  from  (10.12)  - (10.15).  that  for  p sufficiently  small 

(10.16)  RHS  > ^ mAy (t)p.  Ilog  pl,(2|X  < t i K2>  . 

From  (10. 1),  (10.  3)  one  has  again  (10.10)  valid  on  the  interval  being 
considered.  Comparing  (10.10)  and  (10.16)  one  finds  that  (10.  2)  is  satisfied 
on  the  interval  | p |log  p I < t < K^,  with  K2  chosen  as  above,  if  p>0 
is  sufficiently  small.  Note  that  in  this  and  the  previous  part  of  the  proof 
we  have  not  used  the  special  form  of  y(t)  in  any  significant  way. 

Finally,  we  establish  (10.2)  for  any  t > For  any  such  t (10.12) 
is  valid  and  the  estimate  for  given  by  (10 . 13)  holds . Rewrite  ^ in 
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the  form 


llog  m- 


(10.17) 


I.  (w  - <tjt))[  f aft  - s)F  (|<t),|(s))ds] 
1 0 J _ l 


+ (w  - o0(t))[  / aft  - s)F^(  £(t),£  fs))ds]  *11  + \z  ' 


2li 

K 


|log 


Thus  (10.12)  may  be  written  as 

RHS  - 'll  * hz  *lZ*h- 

We  shall  combine  1^  + and  show  that  1^  > |l3l.  From  (10.17),  (10.1), 
(10.13),  and  y nondecreasing,  we  find  easily  that  for  t > K 


I + I > my(K  )|i  flog  (jl  I j a(cr)do- 

LL  L.  £ u 

K2 

. M|mi  {Zy{^2  llog  fx  |2  + M-(g(0)  - b0(0)))  . 

Thus  if  hi  is  sufficiently  small  and  is  not  too  large,  there  exists  a 

constant  6 > 0 such  that  for  t > K 


(10.18) 


!n  + lz  - 6|J-  ^log  K 


I12>|I3I  (t  > k2>  . 


To  show  that 
(10.19) 

it  suffices  to  show  that  for  t > K.,,  filog  (i  | < s < t,  one  has 
(10.  20)  (w  - * HtJFjdft),^))  > (w  - <f0)(s)  |F2(f,(t),i(s))  I . 
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We  write  £(t)  cQ(t)  4 O^w  - cQ)(t),  its)  cQ(s)  + 0^(w  - $ >(sj,  where 

0 < 0^,  0^  < 1 in  the  integrals  I and  in  (10.19).  The  assumption 

2 4 4 

F < C (E  x E and  the  bounds  (10. 5)  insure  that  for  the  range  of  the 
arguments  of  interest  here  the  first  partial  derivatives  of  F^/F’^  do  not 


2M 


exceed  — . Using  these  facts,  as  well  as  F (x,x)  -F  (x  x) 
m 1 2 


for 


every  x,  and  applying  the  mean  value  theorem  to  F /F^,  it  follows 
that  we  will  have  stablished  (10.  20)  (and  also  (10.19)1  if  we  can  prove  that 

(w  ' V(t)  2M 2 

(1°-21)  iw  - cj(si-1+  “1“[  lc0ft)  ' C0fS)i  + lfW  ' C0,(t'  ’ <W  ' 0 S‘  1 


0 


m 


for  t > K and  ^ |log  p | < s < t.  But  for  s,t  so  restricted  it  follows 


2 K 

from  (10.1)  that 


(10.  22) 


(w  - o ) ( t ) - (w  - <jj  )(s)  ...  , . _ 

0 0 y(tj  - ^s;  2 

(w  - <t Q ) ( s ) y(s)  ^ 1 


and  using  (10.  3)  that 
(10.  23)  ^ -^(S}  - exp(K1  / / a(T)drdtr)  - 1 > ^ / / a(T)drdo  . 


t 00 


t 00 


S cr 


S O' 


Thus  using  the  estimate  lo„(t)  - e„(s)  < — / a(o)d  > (obtained  from 

0 0 — m J 

s 

(4.4),  see  also  Remark  5. 1)  and  (10.  23),  in  (10.  21)  it  follows  that  (10.  21) 
holds  if  one  chooses  the  constant  sufficiently  large.  This  proves 


(10.19).  Combining  (10.18)  and  (10.19)  shows  that  for  p sufficiently  small 

0 >0’  Ki 


and  the  constants  6 > 0,  y_  >0.  K,  > 0 chosen  as  above 
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L.f  mathematical  model 


M lten  plastics  commonly  exhibit  large  elastic  recovery:  for  example, 
a filament  'f  a certain  polyethylene  ('Melt  1'  at  150  C , when  elongated 
at  a rate  of  1 cm/sec/cm  from  an  initial  length  of,  say,  1 cm  to  a length 
of  5S  cm  and  then  released,  will  reach  a final  length  of  about  ~>  cm  [1  1 . 
Such  elastic  and  other  rheological  properties  are  of  interest  in  the  process- 
ing of  plastics  and  rubber  and  also  as  examples  of  materials  with  'memory’. 

The  correct  equations  for  describing  the  isothermal  behaviour  of 
a given  molten  plastic  are  not  yet  known.  One  set  which  has  been  used 
with  some  success  [8]  for  'Melt  1'  can  be  expressed  as  follows  in  body 
tensors : 

(Al)  n(P,t)  + py_1(P,t)  -h  + J a(t  - sky  1 ( P , s ids  ; 

_ CO 

(A2)  — det\(P,s)  0 (-*>  < s < t)  : 

dS  “ 

(A3)  ^ • ZL(P,t)  p(£  - H)  ; 

(A 4)  R{Y'P,  s)}  0 (-»  < s < t)  . 

(Al)  is  the  so-called  'rubberlike  liquid'  constitutive  equation  relating 
the  symmetric  contravariant  stress  tensor  iT(P,t)  at  particle  P and 
time  t to  the  values  of  the  reciprocal  y ^ of  the  symmetric  covariant 

The  U.  S.  production  (8  V 10  ' lb)  of  polyethylene  in  ll>73  exceeded  that  of 
any  other  one  polymer. 
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metric  tensor  _V(  P , s ) at  times  in  the  interval  •«<s<t  {see  (A12)  below). 

The  material  properties  are  determined  by  the  nonnegative  constant  r\  and 

the  nonnegative  co  nstants  a, , . . . , a , t , . . . t in  the  'memory  function1 

1’  m 1 m 

m 

(A5i  a ( t ) X a exp(-t/T  ) (a  > 0 t > 0)  . 

— . r r r ~ ’ r 

r 1 

(See  (A26)  and  following  explanation,  below). 

If  the  integral  term  were  omitted,  (Al),  with  the  constant  volume 

condition  (AZ),  would  describe  an  incompressible  Newtonian  liquid  of 

viscosity  q and  could  be  used,  with  the  stress  equation  of  motion  (A3), 

to  derive  the  Navier- Stokes  equations;  there  would  be  no  elastic  recovery 

possible,  a and  H are  contravariant  vectors  describing  acceleration 

and  body  force  per  unit  mass,  respectively;  p denotes  the  density. 

p is  a scalar  function  of  P and  t , of  the  nature  of  a hydrostatic 

pressure,  introduced  in  conjunction  with  the  incompressibility  condition 

(A2).  in  (A3)  denotes  the  covariant  derivative  operator  formed  with 
— t 

jy(P , t) , and  the  dot  denotes  contraction  [8,  p.  193]  . (A4)  expresses 

the  vanishing  of  the  fourth  rank  Riemann-Christoffel  curvature  tensor  R 
constructed  with  y(P,s);  this  expresses  the  fact  that  the  body  manifold 
is  Euclidean  at  each  instant  s and  so  admits  a body  coordinate  system 
that  is  instantaneously  rectangular  Cartesian  at  s [8,  p.  202]  . 

Equations  (Al  - 4)  are  sufficient  in  number  to  determine,  in  principle, 
the  unknown  variables  w , _Y,  and  p when  the  remaining  quantities 
(together  with  suitable  boundary  and  initial  conditions)  are  given. 
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A 


When  referred  to  an  arbitrary  body  coordinate  system  B : fPj  — IRJ, 
(Al)  and  (A2>  yield  equations 


(A6  i , t ) + pfe  ,t  ,t) 


+ J a(t  ' S)ylj(£,,s)ds  , 


(A7 ) — det  y..(4,s)  0 (-s  < s < t)  , 

ds  l]  — ’ 

12  3 

where  4 in  this  and  in  similar  contexts  is  short  for  (4  ,4  ,4  ),  the 

coordinates  of  P in  B;  it11,  \1],  and  y. . denote  components  of  tt. 

i j — ’ 

X \ and  Y,  respectively,  i,j  1,2,3.  For  an  arbitrary  B,  (A3) 
yields  a complicated  set  of  equations.  One  can,  however,  always  choose 
a B that  is  Cartesian  at  the  instant  t at  which  (A3)  applies  (i.e.  , B 

is  such  that  y..(4,t)  is  independent  of  t,  for  i,j  =1,2,3),  and  a 

3 

space  coordinate  system  C : (Q)  — R that  is  rectangular  Cartesian; 

(A3)  then  yields  the  following  three  partial  differential  equations: 

mi  it  i’t^  = p(i¥au)-xi)  ,i  = i,2,3). 

j , k 1 34  34^  9t 

The  motion  of  the  body  is  now  described  by  the  three  equations 


x1  = fx(4 , s; 


(i  = 1,2,3)  , 


where  x1  denote  the  coordinates  in  C of  the  place  0 occupied  by 
particle  T at  time  s.  X*  denotes  the  components  in  C of  the  external 
body  force  per  unit  mass. 

(A 4)  yields  a very  complicated  set  of  nonlinear,  second  order  partial 
differential  equations  in  y_(4,s)  whose  solution  is  of  the  form 
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(A10I 


k 


(i,  j 1,2,3) 


V6'sl 


y £C^.,s) 

^ i 

k 1 0| 


,k 

tit  fc.s) 


for  arbitrary  B and  rectangular  Cartesian  C [ 3 , 11]  . We  may  thus  use 

the  three  functions  f X ( 4 , s ) in  place  of  the  six  functions  y_(4,  s)  ( = -y  (4  , s ) ) 

as  unknowns.  On  using  (A10)  and  the  matrix  equation 


(AID  [ \ij(£ , s)]  [v.At,s)]~l  , 

we  may  express  (A6)  and  (A7)  in  terms  of  tt1-1,  f1,  and  p;  on  substituting 
the  resulting  expressions  for  rr1]  into  (A8),  we  finally  obtain  three  non- 
linear, partial-integro-differential  equations,  which,  with  the  single 
equation  resulting  from  (A7),  yields  a set  of  four  equations  for  the  four 
unknown  functions  f1,  p;  the  independent  variables  are  c,1 , s. 

The  final  equations  are  nonlinear  in  f1,  p although  the  rubberlike 
liquid  constitutive  equation  (Al)  is  linear  in  the  tensors  tt,  y *.  The 
nonlinearity  comes  from  the  constant  volume  condition  (A 2)  and  from  the 
zero -curvature  condition  (A4)  whose  solution  (A10)  is  quadratic  in  the 
unknown  functions  f\  The  nonlinearity  arising  from  the  products  in  the 
left-hand  side  of  (A8)  can  be  removed  trivially  by  choosing  B to  coincide 
with  C at  time  t,  so  that  f^^t)  = 41  and  9fV9£J  = 6. . at  time  t. 

A very  considerable  simplification  of  the  above  equations  is  obtained 
for  flow  histories  which  are  homogeneous  (or  uniform)  under  conditions  in 
which  the  inertial  and  body  force  terms  on  the  right-hand  side  of  (A3)  can 
be  neglected.  Such  histories  are  of  little  or  no  interest  in  classical 
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hydrodynamics  (where  the  constitutive  equations  are  given  by  (Alj  with 
aft)  ; 0)  but  are  of  fundamental  importance  in  polymer  rheology  where 


highly  viscous  molten  plastics  can  be  subjected  to  uniform  elongation  in 
filament  form  or  to  two-way  stretching  in  sheet  form;  results  of  carefully 
controlled  experiments  of  this  type  can  be  used  to  test  the  applicability  of 
constitutive  equations  such  as  (Al).  A flow  history  is  homogeneous  if, 
for  any  two  instants  s,  t,  we  have  _/t_Y(P,s)  0 [8,  p.  247],  from 

which  it  can  be  shown  that  V Y ^(P.s)  0;  since  Y’  commutes  with 

Tf  t 

the  operators  a/9t  and  J aft  - s)...ds,  it  follows  from  (Al)  (taking 
p to  be  independent  of  P,  as  a trial  solution)  that  ( P , t ) 0 (show- 

ing that  the  stress  is  homogeneous)  and  hence  also  that  (A3)  (with  the 
right-hand  side  zero)  is  satisfied.  It  also  follows  from  the  above  defini- 
tion of  a homogeneous  flow  history  that  a body  coordinate  system 
B : P — | exists  that  is  Cartesian  in  every  state,  i.e.  , which  is  such 
that  y..(|,s)  is  independent  of  £ for  all  s [8,  p.  247],  and  hence 
(A4)  is  satisfied.  The  behaviour  in  homogeneous  flow  histories  with  inertial 
and  body  forces  neglected  is  thus  governed  by  (A6)  and  (A7)  with  £ 
absent,  and  there  is  no  longer  a need  to  introduce  a space  coordinate 
system  or  the  functions  f* : one  can  instead  use  y.  or  y1J  as  the 

unknowns,  for  example,  in  the  case  of  problems  involving  the  calculation 
of  free  elastic  recovery:  in  such  problems,  the  flow  history  (and  hence 

v ( s ) ) would  be  specified  throughout  some  interval  -00  < s < t. , say  ; 

'ii  1 
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for  tj  < s < t,  the  stress  would  be  zero,  and  the  elastic  recovery  would 
be  determined  by  solving  the  set  (A6)  (with  ir1]  0)  and  (A7)  for  p(t) 
and  y.ds)  (t^  < s < ti.  These  equations  thus  form  a simultaneous  system 
of  nonlinear  Volterra  integral  equations,  in  which  the  nonlinearity  arises 
from  the  incompressibility  condition. 

In  this  paper,  we  consider  the  particular  case  of  the  above  in  which 
the  specified  flow  history  is  one  of  simple  elongation  (at  constant  volume). 
The  variable  p can  be  eliminated,  and  the  recovery  behaviour  is  then 
governed  by  a single  nonlinear  Volterra  equation,  which  we  now  derive. 

In  any  B,  the  separation  PQP  at  time  t between  any  two 
neighbouring  particles  PQ,  P is  given  by  the  equation 

<A12>  (p0p»t2  - EE 

i j 

where  e,  and  f 6!1  are  the  coordinates  in  B of  P^  and  P. 

For  any  two  times  s,  t,  there  are  three  material  lines  through  any  given 
PQ  that  are  mutually  orthogonal  at  s and  at  t;  in  the  strain  s - t, 
infinitesimal  material  line  elements  tangential  to  these  three  material 
lines  at  PQ  change  in  length  by  factors  X.(PQ,s,t)  which  are  given  by 
the  positive  roots  in  X of  the  equation 

(M3)  d8t{VlJ(|,t)  - - 0 . 

The  factors  X.  are  called  'principal  elongation  ratios ' . A flow  is 
'shear  free'  if  there  exists  a body  coordinate  system  B that  is  always 
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r 


orthogonal,  i.e.,  such  that  y_  0 when  i*j  [8,  p.  81].  Fora 
shear- free  flow,  the  principal  elongation  ratios  are  given  by  the  length 
changes  of  the  coordinate  lines,  and  the  roots  of  (A13)  are  given  by 

(A14)  \Zy..(^,s)  -y . . (t, , t ) (i  1,  Z 3)  . 

l 'll  ’ ’ll 

For  a shear-free  flow  that  is  homogeneous,  B is  always  Cartesian,  i.e.  , 

the  \ are  independent  of  L . A shear-free  flow  is  a 'simple  elongation' 
n 

if  two  principal  elongation  ratios,  X.  and  say,  are  always  equal; 

the  coordinate  lines  are  then  called  directions  of  elongation.  The 
constant  volume  condition  (A7)  reduces  to  the  equation  X \ ^ 1 and 

hence,  for  simple  elongation  at  constant  volume,  we  have 


(A15)  \2  \ = V 


We  now  consider  the  following  problem: 


-- 5 ^ - b 

(A16) 


-t0  < = < » 


(A17 ) 


Zero  stress;  no  flow;  B rectangular  Cartesian;  hence 

ij  n ij  c 

it  0 v, . y c .. 

’ l]  l) 

Homogeneous  simple  elongation  at  constant  volume  and 
constant  rate  k,  i.e., 


dxj_(-t0,  s) 
ds 


“VV”-  VS  S* 


0 < s < t 


Zero  stress;  free  elastic  recovery: 


(A18 


rr1J  0 . 


And  hence  \ are  independent  of  PQ  ■ 
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A 


We  wish  to  calculate  \ (6,s)  for  0 < s < t.  As  a trial  solution, 

ij 

it  is  reasonable  to  suppose  that  the  elastic  recovery  will  involve  a homo- 
geneous simple  elongation  at  constant  volume  with  the  (^-coordinate  lines 
again  as  directions  of  elongation  (or  contraction).  For  convenience,  we 
write 

(A  19 ) y(s)  : x1(-t0»s)  • 

Since  the  entire  flow  history  is  a homogeneous  simple  elongation  at 
constant  volume,  we  have,  from  (All),  (A14),  (A15),  (A16),  and  (A19), 


11..  , -2,  22  . 33  , . . 

y (4,  s)  y (s),  Y (4,s)  = Y (£,s)  = y(s), 


(A  20) 


i] 


(-00  < s < t) 


0 (i  * j ) . 


From  (A16)  and  (A17),  we  have 

S' 

1 


(A  21) 


y(s) 


(-oo  < g < -t  ) 


exp{x(s  + t )}  (-tQ  < s < 0)  . 


y (s ) is  to  be  calculated  for  0 < s < t so  as  to  satisfy  (A6)  and  (A7). 

(A7 ) is  satisfied  by  (A20).  Using  (A20)  and  (A18),  the  six  equations 
(A6)  for  t > 0 reduce  to  the  following  two: 

f 

-2, 


- 2 t 

(A22)  (i  = j = 1)  p(t)y'2(t)  = -t,  dY.  ® 

UL  _ 00 


+ J a(t  - s)y  (s)ds  , 


(A23)  (i  j 2,3)  p(t)y(t)  = + J a(t  - s)y(s)ds  . 

_ 00 


The  unknown  function  p(t)  may  be  eliminated  by  multiplying  (A22)  by  y (t) 
and  then  subtracting  (A23).  The  resulting  equation  may  be  written  in  the  form 


-61- 


(A24)  ' ,Jl  '^dt^  J a(t  - s)F(y(t),y(s))ds  (t  > 0)  , 

- JO 

where 

3 2 

(A25)  p : 3rj , F(y,z)  : (y  /z  ) - z . 

These  are  the  equations  used  in  the  text  above.  Finally,  we  add  some  brief 
remarks  about  the  physical  basis  and  applicability  of  the  rubberlike- 
liquid  constitutive  equation  (Al)  which  has  been  discussed  elsewhere 
(8,  pp.  143,  223-236]  . 

(Al)  has  been  derived  from  two  different  molecular  theories:  the 
'bead-spring'  theory  of  Rouse  and  Zimm  for  very  dilute  solutions  of  deform- 
able long  molecules  in  an  incompressible  Newtonian  solvent  of  viscosity 
rq , and  the  network  theory  of  Green  and  Tobolsky,  Yamamoto,  and  Lodge 
which  is  developed  for  concentrated  polymer  solutions  and  undiluted  or 
molten  polymers.  It  is  curious  that  two  different  molecular  theories 
should  yield  constitutive  equations  of  the  same  form,  but  the  reason  for 
this  is  known:  the  differences  between  the  two  sets  of  equations  at  the 
molecular  level  do  not  survive  the  averaging  process  used  to  go  from  the 
molecular  quantities  to  the  macroscopic  quantities  tt  and  V.  which 
appear  in  the  constitutive  equation  (Al)  [9]  . The  memory  function  constants 

a t are  specified  in  terms  of  three  unknown  constants  by  the  bead- spring 
r r 

theory  but  are  not  specified  by  the  network  theory,  which  also  leaves  q 
unspecified. 
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According  to  the  network  theory,  the  integral  term  in  (Al)  arises  from 
the  thermal  motion  of  a network  composed  of  long,  deformable  polymer 
molecules  temporarily  linked  together  at  a few  points  called  entanglements 
or  temporary  junctions  which  are  assumed  to  be  created  and  lost  at 
constant  rates  which  are  unaffected  by  the  flow  history.  The  concentration 
N(t)dt  of  network  strands  which  were  created  in  the  interval  (0,dti 
and  are  still  in  the  network  at  time  t (a  'strand'  being  that  part  of  a 
polymer  molecule  lying  between  two  consecutive  junctions)  is  given  by  an 
equation  of  the  form 

m 

(A 26)  N(t)  S C exp(-tA  ) 

r 1 r 

where,  for  simplicity,  it  has  been  assumed  that  the  set  of  all  strands 

can  be  sorted  into  m subsets,  labelled  l,2,...,m,  such  that,  in  the 
th. 

r subset,  all  strands  were  created  at  the  same  rate  C (per  millilitre) 
and  have  the  same  probability  1/t  per  second  of  leaving  the  network 
[10]  . The  memory  function  in  (Al)  is  given  by  the  equation  aft)  kTN(t), 
where  k is  Boltzmann's  constant  and  T is  the  absolute  temperature. 
Thus  a - kTC  > 0,  and  (As)  is  proved. 

According  to  the  network  theory,  then,  it  follows  that  aft)  > 0, 
because  there  is  always  a nonzero  concentration  of  strands  of  age  t, 
and  that  a'(t)  < 0 because  strands  of  age  t(>  0)  can  only  be  lost  and 
not  created;  strands  are  created  with  age  0 only.  It  also  follows  from 
(A5)  that 
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(AZ7 ) 


(k  0,1,2,...) 


. . k (k),  _ 

( -1)  a (t)  > 0 ; 

and  a*k*l^t)/a^',it)  is  nondecreasing 


a^'(t)  denotes  the  k^  derivative  of  aft).  (A27)  represents  the 
properties  of  aft)  some  of  whicn  are  used  in  the  present  analysis. 


The  constitutive  equation  obtained  by  putting  g 0 in  (Al)  leads 
to  the  ’reduced'  equation  (1.6)  and  has  been  tested  for  Melt  1 by  comparing 
the  predictions  with  results  of  a series  of  experiments  performed  by 
Meissner  [ 13, 15]  . The  constants  in  (A5),  with  m 5,  were  chosen  to 
fit  stress  growth  data  in  simple  elongation  at  low  rates.  (Al)  (with  g 0) 
then  gave  good  agreement  with  stress  growth  data  in  simple  elongation  at 
higher  rates,  with  elastic  recovery  data  following  elongation  and  following 
shear,  and  with  stress  growth  data  in  shear  flow  | B,  pp.  225-231;  2;  4]  , 
provided  that  the  total  strain  from  rest  was  limited  to  moderate  values; 
at  higher  strains,  there  was  serious  disagreement  between  theory  and  data: 
the  predicted  stresses  and  the  predicted  recoveries  were  greater  than  the 
observed.  The  present  analysis  of  the  elongational  recovery  problem  shows 
that  inclusion  of  the  term  in  g leads  to  a reduction  in  the  predicted 
recovery,  which  is  in  the  right  direction  to  give  better  agreement  with 
experiment.  The  term  in  g has  been  added  to  represent  the  possible 
effect  of  the  presence  of  a viscous  solvent  (in  the  case  of  a concentrated 
polymer  solution)  or  of  low-molecular-weight  polymer  (in  the  case  of  a 
molten  polymer).  The  possible  effects  of  such  a term  are  also  of  some 
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interest  in  connection  with  certain  'fast-strain'  test?  of  the  Gaussian 


network  hypothesis  which  have  recently  been  proposed  as  a possible 
method  of  testing  certain  of  the  network  theory  assumptions  when  separated 
from  the  others  [8,  pp.  231-  236;  12]  . It  is  recognized,  however,  that 
other  modifications  to  (Al)  are  required  if  better  agreement  between  all 
the  predictions  and  data  referred  to  above  is  to  be  obtained  [ 14]  . It 
should,  perhaps,  be  added  that  the  homogeneous  elongation  with  neglect 
of  inertial  and  body  forces  treated  in  the  present  analysis  represents  a 
reasonable  idealization  of  the  conditions  of  Meissner's  elongation 

5 

experiments  on  'Melt  1':  a long  filament  of  high  viscosity  (5  X 10  poise)  was 
floated  on  a bath  of  an  inert  oil,  and  the  homogeneity  of  elongation  was  always 
checked  by  weighing  samples  into  which  the  filament  was  cut  after  elongation; 
the  variation  of  elongation  ratio  along  the  filament  was  about  3%  or  less  [25]  . 

The  present  analysis  depends  on  the  properties 
(A 28 ) F(x,x)  0;  3F(x,z)/3x  > 0;  3F(x,z)/3z  < 0 , 

of  the  function  F(y(t),y(s))  in  the  integrand  of  (1.1).  We  have  been  able 
to  show  that  there  is  reason  to  believe  that  these  properties  would  be 
exhibited  by  the  corresponding  integrands  for  the  same  elastic  recovery 
problem  when  members  of  a much  wider  class  of  possible  constitutive 
equations  than  (Al)  are  considered,  namely,  the  class  of  constitutive 
equations  of  the  form 

^ £ 

(A29 ) Tr(t)  ■ ^(t)  + p(t)6_  -r|  ' v(t)  + / a(t-s)  ^ W^d^^^fSjtlds  , 

-oo  k 
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where  the  common  argument  P has  been  omitted  for  brevity,  Wk  are 
arbitrary  functions  of  the  strain  invariants 

(A30)  IL  11(S)  : _V( t ) ; I - V "V)  : V(s)  , 

b denotes  the  right -covariant  unit  tensor,  and  0 = 0 • 0.  • • • 0 (k 

factors)  where  0 denotes  a right-covariant  strain  tensor  defined  by  the 

equation 

(A3l)  0{s, t)  = V (s)  • V(t)  . 

(A1  is  obtained  from  (A 29)  on  putting  = i,  all  other  - 0,  and 

contracting  from  the  right  with  \-1(t).  (A 29)  has  been  considered  by 
others  [ 1,  5 ] . With  W and  W arbitrary,  no  extra  generality  is 
obtained  by  including  terms  with  k > 2,  because  higher  powers  of  0 
can  be  expressed  in  terms  of  lower  powers  by  a tensor  analogue  of  the 
Cayley -Hamilton  theorem  (remembering  that  the  third  invariant  det  0 is 
equal  to  1 because  of  the  assumed  incompressibility). 

It  is  a straightforward  matter  to  show  that  (A 29)  also  leads  to  an 
equation  of  the  form  (1.1)  with 

(A32)  F(x,  z)  = xG(x/z),  G(a)  - £ W (a)(aZ^  - a ) . 

k * 

Thus  F(x, x)  = 0,  as  required,  and  the  two  remaining  properties  in  (A28) 
will  hold  if,  and  only  if,  G '(a)  > 0 (a  > 1).  If  we  now  put  rj  = 0 and 
consider  the  particular  'step-function"  simple  elongation  history  (at 
constant  volume)  for  which  the  material  is  at  rest  at  zero  stress  up  to  the 


-66- 


instant  0 when  it  is  given  an  instantaneous  simple  elongation  of  magnitude 


a > 1,  it  is  easy  to  obtain  the  following  expression  for  the  tensile  stress 
immediately  after  the  elongation: 

00 

(A33)  (p  - p ) = G(a)  J a(u)du  . 

11  22  0+  -'g 

p denote  Cartesian  components  of  stress  referred  to  a rectangular 
n 

Cartesian  coordinate  system  with  the  1-axis  parallel  to  the  direction  of 
elongation.  Since  it  is  reasonable  to  regard  the  material  in  these 
circumstances  as  (momentarily)  a perfectly  elastic  isotropic  solid  during 
the  instantaneous  elongation,  it  is  reasonable  to  expect  that  the  tensile 
stress  will  increase  monotonically  with  the  elongation  ratio  a.  Since 
the  integral  in  (A33)  is  a positive  constant,  it  then  follows  that 
G'(a)  >0  (a  > 1).  This  completes  the  justification  of  (A28)  for  constitutive 
equations  of  the  form  (A  29). 


-67- 


I 

t 


Appendix  (Bi 

Our  purpose  is  to  collect  the  results  on  integrodifferential  inequalities 
needed  in  the  proof  of  Theorem  8.  These  are  rather  similar  to  classical 
results  of  this  type  for  ordinary  differential  equations  and  Volterra  equations 
developed  e.g.  in  [6]  and  [ 19]  . However,  our  situation  is  sufficiently 
different,  so  that  it  is  simpler  to  give  an  independent  short  exposition 
of  what  is  used,  rather  than  to  try  to  apply  the  known  results. 

In  what  follows  let  D uft)  denote  the  lower  left-hand  Dini 
derivative,  D u(t)  denote  the  upper  left-hand  Dini  derivative,  D u(t) 

-f 

the  lower-right  hand  Dini  derivative,  and  D uft)  the  upper-right  hand 

-f 

Dini  derivative  of  a continuous  function  u.  When  D uft)  ^ D^u(t), 
we  denote  this  common  value  by  u' ft),  the  right-hand  derivative  of  u 
and  similarly  u'  ft)  denotes  the  left-hand  derivative  of  u. 

It  will  be  convenient  to  write  the  initial  value  problem  (1.1),  (1.2) 
in  the  form 


(Bl) 


where 


(B2) 


-py'(t)  / aft  - s) Ffy(t),  y( s))ds  + fg(y)(t) 

0 

y(0)  = gfO)  (p  > 0:  0 < t < °°)  , 


f (y)(t)  = f aft  - s)F(y(t),g(s))ds  . 


Throughout  this  appendix  we  shall  assume 

(B3)  aft)  > 0,  a i/fO,*),  F satisfies  (Hp),  g satisfies  (H  ) . 


1 
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The  basic  result  needed  is  (compare  Theorem  1.  2. 1 of  [ 6]  ). 

Proposition  IB.  Let  the  assumptions  (B3|  be  satisfied.  Let  v,w«  C[[0,y);IR), 
\ > 0 , be  given  functions  satisfying  the  following  properties : 

(B4)  v(0)  < w(0)  , 

r t 

-pD  v(t)  > J a(t  - s)F(v(t),  v(s))  + f (v) ft)  (0  < t < y)  , 

0 a 

( B 5 ) < 

t 

-pD  wft)  < f a(t  - s )F(w(t) , w(s) ) + f (w)(t)  (0  < t < y)  . 

V-  ' 0 a 

Then  v(t)  < w(t)  (0  < t < y). 

Proof.  Define  the  set  Z { t « [ 0,  y)  : w(t)  < v(t)} . If  Proposition  IB  is 
false  the  set  Z * t ; let  ^ Inf  Z.  By  (B4),  t inf  Z > 0 and 

(B6)  v(t^)  w(tj)j  v(t)  < w(t)  (0  < t < t ) . 

Taking  h < 0,  |h  | small,  one  has  v(t^  + h)  < w(t^  + h)  and 

v(t1  + h)  - vltj)  w(t^  + h)  - w(tj ) 

h > h • 

Taking  the  limit  inferior  as  h 0 this  implies 

D_ v(tj)  > D w(tj)  , 

and  therefore, 

-pD  v(t^)  < -pD  w(t^)  (p  > 0)  . 

Applying  this  and  v(t^)  - w(t  ) in  inequalities  (B5)  yields  the  inequality 


(B7)  / a(t^  - s)F(v(tj), v(s)ds  < j a(tj  - s)F(w(tj),w(s))ds  . 
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On  the  other  hand,  the  definition  of  ^ and  (B6),  together  with 
the  assumption  < 0 in  (H  ) implies  that 

F(v(t1) , v(s))  > F(w(tj) , w(s))  (0  < s < t ) , 

so  that,  since  a (t)  > 0, 

*1  h 

<BSl  f a(tA  - s)F(v(t  ),  v(s))ds  > f a (t  - s)F(w(t  ),  w(s))ds  . 

0 0 1 1 

Thus  (B8 ) contradicts  (B7)  and  the  set  Z = 0.  This  proves  Proposition  IB. 

Definition.  We  shall  say  that  w is  an  upper  solution  of  the  initial  value 

problem  (Bl)  on  0 < t < y 11  and  only  if  w e C[  0,  y),  w(0)  > g(0),  w'  (t) 

exists  on  (0,y)  and 

t 

-pw'  (t)  < J a(t  - s)F(w(t),w(s))ds  + f (w)(t)  (0  < t < v)  . 

0 a 

A similar  definition  holds  for  a lower  solution  with  the  inequality  reversed. 
We  remark  that  as  a consequence  of  Theorem  2 the  solution  of  the 

reduced  equation  (1.6)  (recall  that  it  was  proved  in  Theorem  2 that 
<fg(t)  < 0 (0  < t < oo  ))  is  a lower  solution  of  (1. 1)  (or  (Bl) ) on  0 < t < oo  . 

The  main  result  for  the  application  in  Theorem  8 is: 

Proposition  2B.  Let  the  assumptions  (B3)  be  satisfied.  Let  d-  be  the 
solution,  let  w be  an  upper  solution,  and  let  v be  a lower  solution 
of  (Bl)  on  0 < t < y.  Then 

(B9)  v(t)  < cf ( t ) < w(t)  (0  < t < Y)  . 


In  Theorem  8 one  takes  v(t)  cQ(t),  when  <+>0  is  the  solution  of  (1.6), 
and  one  shows  that  w given  by  (10.1)  is  an  upper  solution. 

Proof  of  Proposition  2B.  We  shall  prove  the  second  inequality  in  (B9); 
the  first  is  proved  in  a similar  way.  If  w(0)  > g(0)  the  result  follows 
directly  from  Proposition  IB  with  v replaced  by  4>,  D v by  <£>',  and 
D w by  w1  (t)  (and  one  uses  Lemma  1. 2.  2 of  [ 6] ). 

If  w(0)  = g(0)  one  has  from  (Bl) 

0 

-|x<J»'(0+)  = f (g(0))(0)  / a(-s)F(g(0),  g(s))ds 

a _ ao 

and  from  the  definition  of  upper  solution 

0 

-(J.W1  (0  ) < f (w(0))(0)  J a(-s)F(g(0),g(s))ds  . 

+ 3 -00 

Therefore 

w;(o+)  > <t,'(o+)  , 

and  there  exists  e > 0,  such  that  the  function  m defined  by 
m(t)  = w(t)  - 4>ft)  is  strictly  increasing  on  0 < t < e.  Thus  w(e)  > <fc(ej, 
and  applying  Proposition  IB  as  above  on  the  interval  [ e,y)  shows  that 
for  any  e > 0,  <j;(t)  < w(t)  (e  < t < y).  This  completes  the  proof  of 
Proposition  2B. 
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